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ABSTRACT 


Experimental  and  published  results>  on  the  flow  of  gas 
through  consolidated  porous  media,  have  been  correlated  on  a  log 
log  plot  in  terms  of  the  permeability,  porosity  and  an  inertial 
resistance  coefficient  ^  .  This  plot  indicates  that  an  assumed 
linear  relationship  between  the  three  factors  is  justified  based 
on  the  experimental  data  available.  An  equation  for  this  relation¬ 
ship  is  proposed. 

The  limited  data  obtained,  for  the  flow  of  a  gas  through 
a  porous  media  when  the  temperature  of  the  media  and  of  the 
flowing  gas  was  varied,  indicate  that  the  permeability  is  a 
function  of  the  temperature.  An  early  investigation  reported  in 
the  literature  indicated  that  no  dependence  on  temperature  was 
apparent.  Further  investigation  of  this  effect  needs  to  be  carried 
out. 

Experimental  evidence,  from  the  use  of  large  core  samples, 
points  out  the  need  for  the  proper  design  of  the  equipment  and 
the  use  of  proper  core  mounting  methods.  Significantly  different 
calculated  results  can  be  obtained  on  the  same  porous  core  when 
data  measuring  points  are  incorrectly  located.  An  example 
showing  the  magnitude  of  the  error  caused  by  incorrectly 
located  pressure  measuring  points  is  shown. 
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INTRODUCTION 


A  review  of  the  literature  concerned  with  the  flow  of  a 


gas  through  a  porous  media  was  carried  out.  The  experimental  and 
theoretical  evidence  indicates ,  that  when  the  range  of  validity 
of  Darcy's  law  has  been  exceeded,  an  equation  of  the  form  proposed 
by  Forchheimer  will  describe  the  flow  of  a  gas  through  a  porous 
media.  This  takes  place  when  the  inertial  effects  due  to  the  com¬ 
plicated  nature  of  the  flow  path  become  significant. 


dP 


dL 


A  nomenclature  defining  the  terms  used  in  all  equations  is  shown  at 
the  end  of  this  work.  Equations  similar  to  the  Forchheimer  equation 
have  been  developed  based  on  different  theories.  The  final  form  of 
the  equation  used  in  this  study,  for  the  flow  of  a  perfect  gas 
with  account  taken  of  the  compressibility,  is  the  one  developed  by 
Green  and  Duwez  but  expressed  in  a  slightly  different  form. 


The  terms  in  the  equation  are 

A(P^)  =  difference  of  squared  pressures 

M  =  molecular  weight  of  gas 


g 


conversion  factor 


c 


L 


length  of  porous  sample 


2 


Z 

R 

T 

^  = 
G 

fi  = 


compressibility 
gas  constant 
temperature 
viscosity 
mass  velocity 

coefficient  of  laminar  term  of  flow  equation 
coefficient  of  inertial  resistance  term 
of  flow  equation 


The  factors  o(  and  are  a  measure  of  the  pressure  loss  which 
can  be  attributed  to  the  viscous  resistance  and  inertial 
energy  loss  respectively.  This  equation  is  linear  and  when 
plotted  will  yield  the  values  of  o(  and  /S>  for  the  rock  being 
studied.  The  factor  o(  is  seen  to  be  the  inverse  of  the  permea¬ 
bility  K. 

This  was  the  procedure  followed  during  this  study  in 
which  14  short  core  samples  and  2  long  core  samples  were  used 
to  obtain  experimental  data.  A  similar  procedure  was  used  by 
previous  investigators.  Cornell  and  Katz  obtained  flow  character¬ 
istic  data  on  24  sandstone  samples.  A  log  log  plot  of  the 
flow  characteristics,  K  and  ft  ,  indicated  an  apparent  linear 
relationship  between  them.  They  expressed  this  relationship 
as 


4.11  x  1010 


K  4/ 3 
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A  recorrelation  of  the  data  of  Cornell  and  Katz  by  Janicek  and 

Katz  resulted  in  a  relationship  relating  permeability,  porosity 
and  the  inertial  resistance  coefficient  as 


5.5  x  107 

K5/4  03/ 4 


which  resulted  in  a  set  of  porosity  lines  on  a  log  log  plot  of 
K  against  ^3  .  One  of  the  aims  of  this  investigation  was  to 
attempt  to  justify  or  disprove  the  use  of  such  a  relationship 
with  additional  experimental  data.  An  immediate  observation 
was  that  the  porosity  lines  proposed  by  Janicek  and  Katz  did 
not  appear  to  correlate  the  experimental  data  with  any  consis¬ 
tency.  For  this  reason  the  basis  of  the  correlation  was 
reexamined,  and  an  alternative  correlation  proposed  which  fits 
the  experimental  data  with  much  more  consistency. 

The  permeability,  as  determined  from  a  Klinkenberg 
plot,  is  compared  to  that  determined  from  a  flow  data  plot, 
for  several  of  the  rock  samples  studied.  The  agreement  is 
quite  good,  but  there  appears  to  be  justification  for  obtaining 
flow  data  plots  at  different  constant  mean  pressures.  Such 
plots  would  ensure  that  what  is  assumed  to  be  the  straight  line 
flow  region  of  the  flow  data  plots  is  not  being  affected  by 
gas  slippage. 

The  dependence  of  the  permeability  on  temperature  was 
examined  by  Johnston  and  Taliaferro  for  the  flow  of  gas  through 
consolidated  porous  cores.  Their  findings  indicated  that  no 
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such  dependence  was  apparent.  Greenberg  and  Weger,  when  inves¬ 
tigating  the  possible  dependence  of  and  K  on  temperature  for 
gas  flow  through  porous  sintered  metals,  found  that  /3  was  un¬ 
affected  but  that  with  increased  temperature  the  permeability 
decreased.  It  was  felt  that  increased  temperature  might  have 
an  effect  on  the  permeability  of  a  consolidated  porous  rock  and 
for  this  reason  a  series  of  experimental  runs  were  carried  out 
on  one  of  the  small  cardium  sandstone  samples.  The  results 
appear  to  agree  with  the  findings  of  Greenberg  and  Weger. 
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THEORY 


The  flow  of  gas  through  porous  materials  has  been 

extensively  investigated  both  theoretically  and  experimentally. 

Molecular  streaming,  viscous  flow  and  non  Darcy  or  "turbulent" 

flow  are  the  regions  of  flow  which  have  been  considered 

separately.  An  excellent  summarization  of  the  significant 

theoretical  and  experimental  work  carried  out  in  the  three  flow 

( 27 ) 

regions  has  been  compiled  by  Scheidegger' 

Deviations  from  Darcy's  law  have  been  observed  and  studied 

at  low  flow  rates.  As  the  radius  of  the  flow  channels  becomes 

smaller  the  viscous  flow  of  the  gas  will  undergo  a  transition 

to  slip  flow  and  thence  to  molecular  streaming.  This  flow 

regime  has  been  successfully  described  by  Grunberg  and  Nissan^^, 
(21)  (7) 

Klinkenberg  and  Carman'  '  .  When  the  mean  free  path  of  the 

gas  molecules  approaches  the  dimensions  of  the  flow  path,  the 
rate  of  flow  exceeds  that  calculated  by  Darcy's  law.  This 
increase  is  attributed  to  slippage  at  the  walls  of  the  pores. 

At  very  low  pressures  the  mean  free  path  will  be  greater  than 
the  dimensions  of  the  flow  path  and  viscosity  will  no  longer 
be  a  factor  in  the  flow  process.  The  gas  molecules  will  collide 
only  with  the  walls  of  the  flow  path  and  not  with  each  other; 
leading  to  a  form  of  diffusion  flow. 

As  flow  rates  are  increased,  a  deviation  from  viscous 
flow  takes  place.  Opinion  varies  as  to  the  mechanism  of  the 
energy  loss  causing  the  deviation  and  as  a  consequence  the 
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terms  "turbulent",  "quasi  turbulent",  "non  Darcy"  and  "inertial" 
have  been  applied  to  the  flow  region.  It  has  been  shown  that 
the  dependence  of  the  pressure  drop  on  the  flow  rate  becomes 
non  linear  as  soon  as  the  pressure  loss  due  to  inertia  becomes 
significant ) (18) (27) ^  in  curVed  tubes  this  is  long  before 
the  onset  of  turbulence.  Numerous  authors  have  proposed 
theories  based  on  their  concept  of  the  energy  loss  mechanism 

in  porous  media(^)  (^)  (9)  (17)  (24)  . 

( 17 ) 

Hubbert'  '  proposed  an  "inertial  resistance"  where  the 
distortion  of  the  streamlines  due  to  the  inertial  forces 
increases  the  velocity  gradient,  which  increases  the  frictional 
forces  over  what  they  would  have  been  had  no  distortion  occurred 
Green  and  Duwez^^^  support  the  theory  that  distortion  of  the 
flow  lines  and  the  possible  formation  of  stationary  vortices 
can  be  expected.  Tek,  Coates  and  Katz  ^2)  state  that  the  abrupt 
transition  from  streamline  flow  to  turbulent  eddies  has  been 
observed  by  one  of  the  authors  in  media  with  large  pores  but 
conclude  that  the  irreversibilities  occur  due  to  extra  fluid 
motion  consuming  energy,  probably  in  the  absence  of  turbulent 
eddies.  Turbulence  can  be  defined  as  the  random  variation  with 
time  of  the  velocity,  at  any  point  in  the  flow  stream.  In 
terms  of  this  definition  the  general  opinion  is  that  true 
turbulence  is  unlikely  in  fine  grained  porous  material.  The 
initial  deviation  from  a  linear  relationship  undoubtedly  takes 
place  in  the  laminar  flow  region  when  inertial  effects  become 
appreciable  due  to  the  distribution  of  pore  sizes  with  their 


-  7  - 

rapidly  changing  directions  and  cross  sections.  The  term  "non 
Darcy"  will  be  used  to  describe  the  condition  of  gas  flow  where 
the  difference  of  squares  of  pressures  are  more  than  propor¬ 
tional  to  increases  in  flow  rates. 

Because  of  the  assumption  of  an  analogy  between  the  flow 

of  fluids  in  tubes  and  in  porous  media,  much  effort  has  been 

*  *  t  » 

directed  toward  finding  the  universal  Reynolds  number  at  which 
a  deviation  from  Darcy's  law  will  occur  in  all  porous  media. 
Since  the  Reynolds  number  and  friction  factor  are  both  dimen¬ 
sionless  groups  one  must  be  a  numerical  function  of  the  other. 
However,  the  usual  representation  of  the  flow  system  by  a 
Reynolds  number  -  friction  factor  correlation  depends  on  the 
definition  of  a  characteristic  length  for  the  media,  which  for 
consolidated  media  is  difficult  to  establish.  Such  correla¬ 
tions  have  been  achieved  for  special  systems  ^  ^  and 

are  useful  when  dealing  with  such  systems.  However  there  does 
not  appear  to  be  any  justification  for  assuming  flow  to  be 
analogous  when  the  Reynolds  number  for  different  systems,  are 
the  same.  If  the  porosity,  tortuosity  and  cross  section  were 
the  same,  different  results  would  still  be  obtained  if  the 
curvature  of  the  flow  paths  differed. 

Experimental  evidence  has  shown  that  rather  than  an 
empirical  correlation  between  the  friction  factor  and  Reynolds 
number  an  equation  of  the  type  suggested  by  Forchheimer,  will 

"  SUL  =  Ml  +  bi^2 
dL  K 

correlate  the  flow  data  results  for  one  dimensional  flow. 
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A  nomenclature  listing  and  defining  the  terms  used  in  all 

equations  is  shown  at  the  end  of  this  work.  The  results  of 

attempts  to  develop  an  equation  to  explain  non  Darcy  flow  show 

that  an  equation  of  this  type  should  describe  the  flow  data 

when  inertial  effects  are  significant. 

( 19) 

Johnston  and  Taliaferro  found  the  empirical  equation 


(2) 


would  express  the  resultant  gas  flow  data  curve  for  both 

consolidated  and  unconsolidated  porous  samples.  Green  and  Duwez ' 

developed  a  similar  equation  from  the  dimensional  considerations 

( 24) 

proposed  by  Muskat  as 

-  =  o<  (3) 


where  o(  and  p  define  the  viscous  and  inertial  resistance  of 
the  porous  system.  These  two  factors,  which  have  the 
dimensions  of  length,  were  considered  independent  of  the 
mechanical  properties  of  the  fluid  but  characterize  the  structure 
of  the  porous  system.  The  factor  o(  is  the  inverse  of  the  per¬ 
meability  factor  K. 

The  theories  based  on  the  Kozeny  theory  of  capillaric 
models  are  restrictive  in  nature  and  were  not  normally  applied 
to  consolidated  media.  Ergun  and  Orning^"^  developed  a 
quadratic  equation  describing  flow  through  unconsolidated  media 
based  on  this  theory.  Their  equation  was  developed  from  the 
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Reynolds  theory  of  resistance  to  flow  and  the  Kozeny  assumption 
that  the  media  can  be  considered  analogous  to  a  bundle  of 
parallel  capillaries  with  a  common  hydraulic  radius.  °(  and  (2 
were  used  to  indicate  the  part  of  the  pressure  drop  which  was 
due  to  viscous  and  inertial  energy  loss  and  they  are  necessarily 
different  for  each  media.  Just  as  it  has  not  been  possible 
to  determine  the  Reynolds  number  on  a  universal  basis,  so  are 
o(  and  fi  dependent  on  the  structure  of  the  porous  media. 


A  similar  approach  was  used  by  Cornell  and  Katz^^  to 


obtain  a  quadratic  equation  of  flow  for  consolidated  media. 

A  Kozeny  type  equation  was  written  in  terms  of  the  hydraulic 
radius  and  the  tortuosity,  as  measured  electrically.  The 
problem  of  determining  their  factors  k-^  and  "k.2  on  a  universal 
basis  remains  and  they  have  to  be  determined  for  each  rock 
sample.  The  equation  they  developed  is 


Equation  (3)  proposed  by  Green  and  Duwez^^^  can  be  rewritten 


for  the  flow  of  gas  by  using  the  perfect  gas  law  with  a 
compressibility  coefficient  and  assuming  that  under  steady 
state  conditions  the  mass  velocity  of  the  gas  p  v  - 
constant.  The  resulting  equation  is 


2 


G  is  a 
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When  the  integration  is  carried  out  the  final  equation  becomes 


(5) 


This  equation  when  plotted  yields  a  straight  line  with  an 
intercept  of  c(  and  a  slope  of  ^  .  When  compared  to  the  Cornell 
and  Katz  equation  it  can  be  seen  that 


o(  =  32F 


Equation  (5)  was  the  one  used  by  Cornell  and  Katz  and  by 

( 15) 

Greenberg  and  Whger v  '  in  analysing  their  flow  data.  It  is 


the  equation  that  was  used  in  this  study  to  analyse  similar 
flow  data.  Appendix  II  contains  an  example  calculation,  showing 
the  calculations  necessary  for  each  flow  rate  data  point.  All 
of  the  equations,  with  their  constants  determined  for  the  units 
of  measurement  used  in  this  study,  are  included  with  the 
example  calculation  in  Appendix  II.  The  data,  on  the  physical 
properties  of  gases  and  liquids,  necessary  for  the  calculations 
are  shown  in  Appendix  III. 

By  definition  the  Reynolds  number  of  a  flow  system 
represents  the  ratio  of  the  inertia  forces  to  the  viscous 
forces^^  .  From  equation  (5)  this  can  be  seen  to  be 


(8) 
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Similarly  the  friction  factor  is  defined  by  the  ratio  of  the 
total  pressure  loss  to  the  inertia  energy  loss,  which  from  the 
same  equation  is 


f  =  A(P2)  M  9r. 

ft  L  (ZRT)  G2 

Equation  (5)  can  then  be  written  as 


(9) 


f 


2 


(10) 


For  a  consolidated  system  it  is  impossible  at  present  to 
describe  the  geometry  of  the  system  in  detail.  Factors  such 
as  porosity,  tortuosity,  pore  size  distribution,  flow  path 
curvature,  particle  size  and  roughness  are  just  some  of  the 
factors  which  govern  the  values  of  c(  and  p  .  Green  and  Duwez 
suggest  that  the  friction  factor  correlation  using  equations 
(8)  and  (9),  which  requires  two  characteristic  length  parameters, 
is  the  only  one  that  can  be  successfully  applied  to  consolidated 
materials . 

c(  and  ft  are  usually  assumed  to  be  constant  for  specific 
porous  systems.  The  validity  of  this  assumption  was  investi¬ 
gated  by  Greenberg  and  Weger  for  gas  flow  through  porous  sintered 
metals.  There  was  no  apparent  effect  on  either  coefficient  due 
to  pressure  but  with  increased  temperature  the  value  of  o{  was 
found  to  increase.  This  has  not  been  investigated  to  any  extent 
for  porous  rock  systems.  While  a  temperature  effect  would 
not  be  significant  in  an  essentially  isothermal  reservoir,  the 
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usual  practice  is  to  determine  the  permeability  at  room 
temperature  without  taking  account  of  the  reservoir  tempera¬ 
ture  . 


EXPERIMENTAL  EQUIPMENT 


As  part  of  the  experimental  investigation  it  was  initially 
decided  that  large  core  samples  would  be  mounted  in  lucite  and 
in  steel  cylinder?,  in  which  pressure  taps  could  be  fitted  for 
the  determination  of  pressure  drop  along  the  length  of  the  core 
under  conditions  of  high  and  low  flow  rates.  It  was  considered 
possible  that  the  point  of  deviation  from  Darcy's  equation  could 
be  observed  along  the  length  of  the  core  under  high  flow  rates. 

The  mounting  of  the  cores  proved  difficult.  An  epoxy 
resin  mounting  material,  Shell  epon  828  with  a  filler  material 

/  1  r  \ 

and  curing  agent,  similar  to  that  used  by  Greenburg  and  Weger ' 
was  employed.  This  material  was  generally  unsuccessful  for 
cores  of  greater  than  one  inch  diameter.  Its  apparent 
shortcomings  were  that  it  contracted  on  cooling  leaving  a 
void  space  between  the  solidified  material  and  the  casing. 

Further  additions  of  material  did  not  fill  these  voids  completely 
leaving  permanent  spaces  next  to  the  casing  wall.  Examination 
also  showed  that  considerable  air  was  entrapped  in  the  material 
in  the  form  of  bubbles.  A  further  source  of  difficulty  was  in 
attempting  to  effect  a  seal  between  the  end  plates  and  the 
mounted  core.  The  best  method  appeared  to  be  "0"  rings,  when 
the  end  faces  provided  a  good  seating  surface  and  did  not  have 
large  air  bubbles  exposed. 

Under  high  pressure  the  resin  material  permanently  deformed 
and  cracks  invariably  developed.  Often  these  cracks  were  propa¬ 
gated  into  the  porous  core.  The  holes  drilled  for  the  pressure 
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taps  were  a  point  of  weakness  which  together  with  the  void 
spaces  next  to  the  casing  wall  were  considered  suspect  as 
possible  means  of  pressure  communication.  Flow  data  and 
results  of  such  a  mounting  procedure  are  discussed  in  a  later 
section. 

The  use  of  a  non  shrink  concrete,  used  commercially  as 
a  grouting  material,  was  investigated  as  a  core  mounting 
material.  Any  size  of  core  can  be  mounted  at  one  pouring  but 
extreme  care  is  necessary  to  remove  all  of  the  entrained  air. 
When  the  cement  mixture  sets  the  permeability  was  found  to  be 
less  than  0.1  millidarcies  and  decreased  with  time.  Prior  to 
mounting  cores  in  this  cement  they  were  given  a  thin  coat  of 
the  previously  mentioned  epoxy  resin  material. 

The  major  experimental  data  were  obtained  on  small  core 
samples  of  consolidated  porous  media.  Various  methods  of 
holding  such  samples  are  described  in  the  literature ^  (15) (18) . 
A  modification  of  the  method  described  by  Johnston  and  Talia- 
ferrov  u/ ,was  used  and  is  shown  in  Figure  1.  The  holder  was 
machined  as  shown  from  solid  stock.  The  rubber  stoppers  used 
were  #26  extra  long,  drilled  to  one  inch  diameter,  which  fitted 
in  the  tapered  section  of  the  holder.  Pressures  up  to  1900 
psia  were  imposed  across  the  short  core  samples  without 
difficulty. 

Figure  2  is  a  composite  diagram  of  the  experimental 
apparatus  used  for  obtaining  data  on  both  long  and  short  core 
samples.  Nitrogen  cylinders  were  manifolded  through  a  pressure 
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regulator  to  give  a  sufficient  source  of  gas  for  the  flow 
measurements.  Bourdon  tube  pressure  gauges  were  used  for 
pressure  measurements  greater  than  two  atmospheres  gauge 
pressure.  Below  that  pressure,  water  or  mercury  manometers  were 
used  for  the  measurement  of  pressures.  Temperature  measurement 
of  the  upstream  and  downstream  flowing  temperatures  were 
made  using  iron  -  constantan  thermocouples.  Gas  flow  measure¬ 
ments  were  made  at  low  flow  rates  using  burettes,  from  1  to  25 
cfm  using  a  wet  test  meter  and  above  25  cfm  using  a  diaphragm 
meter.  The  meters  were  calibrated  using  a  gas  prover. 

The  pressure  drop  profile  was  obtained  on  long  core 
samples  using  calibrated  transducers  which  were  recorded  on  a 
continuously  recording  visicorder  plot  as  pressure  drop  per 
section  in  psi.  A  constant  voltage  source  and  digital  voltmeter 
were  used  to  ensure  a  constant  voltage  input  to  the  transducers. 

A  possible  defect  in  the  core  holding  apparatus  used, 
was  that  gas  might  bypass  the  core  in  some  way  during  the  experi¬ 
mental  runs,  either  between  the  rubber  stopper  and  the  holder 
wall  or  between  the  core  and  the  rubber  stopper.  To  help 
prevent  the  gas  passing  between  the  stopper  and  the  holder 
wall  the  stopper  was  greased  before  being  placed  in  the  holder. 

In  an  attempt  to  determine  if  gas  could  bypass  the  core  a 
machined  lucite  plug  was  placed  in  the  holder.  When  a  pressure 
of  60  psig  was  applied  to  the  upstream  section  of  the  core 
holder  no  flow  was  evident  immediately,  then  a  very  slow  seepage 
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was  apparent  which  built  up  to  a  static  pressure  of  4  inches 
of  water  in  three  minutes.  The  flow  rate  was  in  effect  negli¬ 
gible.  It  is  felt  that  this  seepage  may  have  been  due  to  the 
surface  of  the  lucite  cylinder,  which  is  extremely  smooth, 
making  it  an  extremely  poor  gripping  surface.  The  rock  samples 
however  offered  a  rougher  surface  allowing  a  firm  contact 
between  the  rubber  and  the  rock  material.  Care  was  taken  that 
the  follower  in  the  core  holder  was  screwed  down  as  tightly 
as  possible  for  all  of  the  experimental  runs.  As  far  as  could 
be  determined  there  would  not  have  been  any  significant  error 
in  the  flow  data  due  to  the  core  being  bypassed  by  the  gas 
flowing  through  the  core  holder. 


EXPERIMENTAL  PROCEDURE 


Sample  Preparation 

Small  cylindrical  samples  of  the  porous  media  investi¬ 
gated  were  cut  from  bulk  samples  using  a  one  inch  diamond 
drill.  All  samples  were  cut  parallel  to  the  bedding  planes. 

They  were  then  cleaned  of  any  loose  surface  material  and  ex¬ 
tracted  in  Soxhlet  extractors  using  benzene.  The  extraction  was 
continued  for  at  least  12  hours  after  which  time  the  samples 
were  left  to  soak  in  benzene  for  a  further  period  of  12  hours. 

If  after  this  time  the  benzene  was  found  to  be  clear,  the 
sample  was  dried  in  a  constant  temperature  oven  for  three  hours 
at  200°F.  If  the  benzene  was  found  to  be  oil  stained  after 
the  soaking  period,  the  extraction  was  continued  until  the 
benz  ene  was  found  to  be  clear  after  soaking  for  12  hours. 

After  drying,  the  samples  were  allowed  to  cool,  and  kept  stored, 
in  a  dessicator  until  used. 

The  samples  were  cut  to  approximately  ±h  inches  in 
length,  with  the  exposed  faces  produced  by  fracturing.  Figure 
3  shows  the  details  of  the  apparatus  used  to  break  off  the 
ends  of  the  sample.  The  sample  was  held  tightly  in  the  apparatus 
while  pressure  was  applied  to  the  knife  blades.  The  fracturing 
can  be  caused  in  two  ways,  by  sharp  blows  to  the  knife  blades, 
and  by  applying  the  pressure  slowly  to  the  blades  using  a  bench 
vise.  Both  methods  gave  satisfactory  results.  The  softer 
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sandstones  gave  the  best  results  when  the  pressure  was  applied 
slowly  while  the  dense  limestones  and  well  consolidated  sands 
were  successfully  fractured  using  both  techniques. 

Porosity  Measurement 

The  bulk  volume  of  the  samples  was  determined  using  a 
mercury  pycnometer.  The  pycnometer  was  first  calibrated  using 

lucite  plugs  which  had  been  machined  to  the  required  dimensions. 

(2) 

A  modified  Oilwell  Research  Porosimeter  using 
Helium  gas  was  used  to  determine  the  matrix  volume  of  the 
samples.  This  apparatus  is  shown  in  figure  4.  Machined  steel 
blanks  were  used  to  calibrate  the  porosimeter. 

The  average  diameter  of  the  samples  was  determined  by 
averaging  five  different  diameter  measurements  along  the  length 
of  the  sample  as  measured  using  a  micrometer.  The  length  of 
each  sample  was  computed  from  the  data  on  its  volume  and 
calculated  cross  sectional  area. 

The  bulk  volume  as  determined  from  the  mercury 
pycnometer  measurements  and  the  matrix  volume  as  determined 
by  the  Oilwell  Research  porosimeter  were  used  to  determine 
the  porosity  of  the  small  core  samples. 

In  the  case  of  the  long  core  samples  the  pore  volume 
was  measured  by  evacuating  the  weighed  cores,  saturating 
them  with  water  and  reweighing.  The  bulk  volumes  were 
determined  by  measurement  of  their  length  and  cross  sectional  areas. 
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Table  1  is  a  tabulation  of  the  physical  properties  of  the 
small  core  samples. 


Table  1 

Small  Core  Physical  Properties 


Core 

Formation 

Rock  Type 

Core 

Diameter 
(  cm)  

Bulk  Volume 
( Pycnometer) 
(cc) 

Equivalent 

Length 

(cm) 

1 

Cardium 

Sandstone 

2.532 

18.751 

3.723 

2 

Cardium 

Sandstone 

2.535 

18.886 

3.742 

3 

Cardium 

Sandstone 

2.535 

19.094 

3.783 

4 

Cardium 

Sandstone 

2.532 

18.690 

3.711 

5 

Cardium 

Sandstone 

2.530 

18.664 

3.713 

6 

Sturgeon  Lake 

Limestone 

2.522 

19.571 

3.917 

7 

Sturgeon  Lake 

Limestone 

2.531 

19.150 

3.806 

8 

Berea 

Sandstone 

2.529 

18.806 

3.745 

9 

Indiana 

Limestone 

2.529 

18.208 

3.626 

10 

Blairmore 

Basal  Quartz 

2.579 

19.165 

3.670 

11 

Blairmore 

Basal  Quartz 

2.579 

17.793 

3.405 

12 

Viking 

Sandstone 

2.579 

18.757 

3.590 

13 

Viking 

Sandstone 

2.580 

19 . 47  2 

3.726 

14 

Belly  River 

Sandstone 

2.579 

19.629 

3  o  7  58 
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Table  1  -  continued 


Core 

Pore 

Volume 

(cc) 

Porosity 

(%) 

Length 

Inches 

Area 

Inches 

1 

3.656 

19.50 

1.466 

.781 

2 

.446 

2.36 

1.473 

.7  82 

3 

3.593 

18.82 

1.490 

.782 

4 

3.945 

21.11 

1.461 

.781 

5 

3.202 

17.16 

1.462 

.779 

6 

2.848 

14.  55 

1.542 

.774 

7 

2.605 

13.60 

1.498 

.7  80 

8 

3.646 

19.39 

1.474 

.778 

9 

3.543 

19.46 

1.428 

.778 

10 

3.825 

19.96 

1.445 

.810 

11 

3.628 

20.39 

1.341 

.810 

12 

1.842 

9.82 

1.414 

.810 

13 

2.002 

10.28 

1.467 

.810 

14 

3.701 

18,85 

1.480 

.810 

Core  diameter  -  average  of  5  measurements 


Bulk  volume  -  average 
Equivalent  length  = 


of  3  measurements 

Bulk  Volume 
Area 

Pore  Volume 


Porosity- 


Bulk  Volume 
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Flow  Measurement 

Flow  data  were  obtained  on  14  small  core  samples  and  2 
large  core  samples  using  the  experimental  setup  shown  in 
figure  2.  The  small  cores  were  placed  in  the  core  holder  as 
shown  in  figure  1.  The  gas  supply  was  connected  to  the  upstream 
end  of  the  holder  and  the  flow  measuring  equipment  to  the 
downstream  end.  Connections  were  made  for  the  measurement  of 
upstream  and  downstream  pressure.  The  thermocouples  used  for 
the  determination  of  upstream  and  downstream  flowing  temperatures 
were  connected  to  a  Leeds  and  Northrup  potentiometer.  All 
connections  were  tested  under  pressure  with  a  soap  solution, 
particular  care  being  taken  to  ensure  that  all  downstream 
connections  and  pressure  measuring  connections  were  gas  tight. 

Two  flow  measurements  were  taken  at  each  data  point  to 
ensure  that  stable  conditions  of  flow  had  been  established. 

A  period  of  5  to  10  minutes  was  allowed  to  elapse  between  the 
two  measurements.  Any  significant  change  in  upstream  pressure 
meant  a  further  period  of  flow  was  required  to  attain  stable 
conditions.  Only  at  highflow  rates  was  any  difficulty 
experienced  in  attaining  stable  conditions.  This  could  have  been 
minimized  by  having  a  larger  gas  supply  or  a  method  of  main¬ 
taining  the  pressure  of  the  gas  supply  at  a  constant  value. 

No  attempt  was  made  to  maintain  a  constant  upstream 
and  downstream  flowing  gas  temperature  for  the  majority  of  the 
runs.  The  gas  inlet  line  was  \  inch  copper  tube  which  between 
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the  regulator  and  the  inlet  connection  had  a  coil  of 
several  feet  immersed  in  a  water  bath  to  bring  the  temperature 
of  the  gas  up  to  the  room  temperature  after  the  cooling  caused 
by  expansion  through  the  regulator. 

Data  on  one  of  the  core  samples  were  obtained  at  three 
different  temperatures.  A  constant  temperature  water  bath  and 
heating  tapes  were  used  to  stabilize  the  core  temperature  and 
the  flowing  gas  temperature  at  the  desired  value.  Extreme 
difficulty  v/as  experienced  in  stabilizing  the  temperature  at 
the  low  flow  rates.  It  v/as  necessary  to  allow  a  lengthy 
period  of  flow  after  every  minor  heating  tape  adjustment,  to 
ensure  that  a  stable  temperature  condition  was  attained. 

All  measured  flow  data,  and  calculated  flow  characteristic 
data,  are  tabulated  in  Appendix  I  and  are  not  included  in  the 
body  of  this  report.  Tables  referred  to  by  number,  in  the 
discussion  following,  are  located  in  Appendix  I. 
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DISCUSSION  OF  EXPERIMENTAL  RESULTS 

Flow  tests  were  carried  out  initially  on  two  large 
cores.  An  Indiana  limestone  core  was  mounted  in  a  steel 
casing  using  an  epoxy  resin  mounting  material.  Attempts  to 
determine  the  porosity  of  this  core  using  saturation  techniques 
were  unsuccessful.  The  porosity  calculated  using  a  representa¬ 
tive  small  core  sample  was  assumed  correct  due  to  the  homo¬ 
geneity  of  the  Indiana  limestone.  The  results  of  the  flow 
tests  on  this  core  are  recorded  in  table  5  in  Appendix  I  and 
are  shown  on  figure  24.  These  were  flow  tests  carried  out 
prior  to  drilling  pressure  taps  along  the  length  of  the  core. 
The  plotted  flow  data  would  appear  to  indicate  that  the 
straight  line  predicted  by  equation  (5)  had  been  established. 

A  comparison  of  figure  24  with  figure  21,  which  shows  results 
obtained  from  flow  tests  using  a  small  Indiana  limestone  core 
sample,  would  indicate  that  the  flow  data  were  probably  still 
in  the  flow  region  affected  by  molecular  streaming.  This  would 
result  in  a  permeability  value  which  would  be  too  high  and  a 
^0  value  which  would  be  too  high.  While  the  small  core  sample 
may  well  be  of  lower  permeability  than  the  long  core  sample, 
it  would  appear  necessary  to  take  flow  data  over  a  more 
extended  range  of  flow  rates  before  the  permeability  determined 
on  the  long  core  sample  could  be  assumed  to  be  correct.  No 
further  data  was  taken  on  this  core  because  of  the  trouble 
experienced  with  epoxy  resin  mounted  cores. 
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FIGURE  21  FLOW  DATA  PLOT  TO  OBTAIN  <*  ANO 


jj  (Feet*1) 


FIGURE  24  FLOW  DATA  PLOT  TO  OBTAIN  *  AND  ft 


An  Alundum  core  was  mounted  in  a  lucite  cylinder  using 
the  epoxy  resin  mounting  method.  The  main  difficulty/ experienced 
With  this  core,  was  obtaining  a  seal  between  the  end  plates  and 
tHe  mounted  core.  "0"  rings  proved  successful  at  the  pressures 
used.  The  maximum  pressure  applied  at  the  upstream  end  of  the 
core  was  approximately  240  psig.  At  this  and  lower  pressures 
the  "O"  rings  maintained  an  effective  seal  and  no  visible 
cracking  within  the  epoxy  resin  occurred. 

Initial  flow  data  wore  obtaihed  for  a  small  range  of  flow 
rates.  That  data  is  recorded  in  table  4(A) .  After  drilling 
pressure  taps  along  the  iong'h  of  the  core  further  data  runs 
were  carried  out,  these  are  recorded  in  table  4(B) .  Until 
this  time  the  upstream  and  downstream  pressure  taps  were  located 
in  the  flow  line  just  upstream  and  downstream  of  the  end 
plates.  Using  the  pressure  taps  and  recording  apparatus  shown 
in  figure  2  a  record  of  the  pressure  drop  across  each  section 
of  the  core  was  obtained.  Tf  :so  data  were  then  converted  to 
absolute  pressure  at  each  point  and  the  resultant  pressure 
squared  tabulated  in  table  2.  Figure  5  is  a  plot  of  the 
pressure  squared  along  the  length  of  the  core.  The  increasing 
outlet  pressure  squared,  and  the  increasing  slope  of  the 
lines,  indicates  measurements  taken  at  increasing  flow  rates. 

The  decreasing  slope  li.nes,  with  the  inlet  pressure  squared 
remaining  essentially  constant,  are  for  runs  where  an  in¬ 
creasing  back  pressure  was  applied  to  the  core  resulting  in  lower 
flow  rates.  The  flow  rates  for  the  plotted  runs  are  tabulated 


in  table  2. 
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Figure  5  shows  that  a  pressure  drop  in  the  two  end 
sections,  greater  than  would  be  predicted,  was  experienced. 
Because  of  this  result  the  pressure  measuring  points  were  changed 
to  the  end  plates  and  further  flow  data  were  obtained.  The 
resultant  plots  of  the  pressure  squared  along  the  length  of 
the  core  are  shown  in  figure  5  and  figure  6.  Two  separate 
plots  were  necessary  because  of  the  range  of  flow  rates 
covered.  Figure  5  shows  the  low  flow  rate  data  while  Figure 
6  shows  the  high  flow  rate  data.  Straight  lines  resulted  from 
the  plot  of  these  data  which  is  in  accordance  with  the  predicted 
result  for  Darcy  flow.  Very  slight  deviation  from  the  straight 
line  was  observed  at  the  higher  rates  of  flow,  shown  on  Figure 
6,  but  it  was  not  sufficiently  pronounced  to  indicate  a 
definite  deviation  from  Darcy  flow. 

The  relocation  of  the  pressure  taps  indicated  that  a 
straight  line  relationship, of  pressure  squared  along  the  length 
of  the  core,  existed  for  the  range  of  flow  rates  covered.  The 
flow  rates,  when  the  pressure  taps  were  located  in  the  flow 
lines,  were  low  enough  to  have  given  a  similar  linear  result. 

For  this  reason  the  plots  of  those  data,  shown  on  figure  5,  were 
extrapolated  to  give  predicted  inlet  and  outlet  pressures 
squared.  From  these  extrapolated  pressures  squared,  the 
predicted  pressures  at  the  inlet  and  outlet  faces  of  the  core 
were  calculated.  Table  2  shows  the  error  in  the  measured 
pressure  drop,  in  the  two  end  sections,  due  to  the  measurement 
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of  the  inlet  and  outlet  pressure  in  the  flow  line.  The 
error  due  to  the  inclusion  of  the  short  lengths  of  flow  line 
is  seen  to  be  of  significant  magnitude. 

Also  shown  in  figure  6  is  the  result  that  was  obtained 
when  the  length  of  line,  from  the  end  plate  outlet,  to  the 
gas  flow  measuring  apparatus  was  shortened.  The  reduction  in 
back  pressure  allowed  a  higher  flow  rate  to  be  acheived  for  the 
same  upstream  or  inlet  pressure. 

The  plotted  flow  data  for  the  Alundum  core  is  shown  in 
figure  7.  Curve  A  represents  the  condition  where  the  pressure 
taps  we re  in  the  flow  line,  for  curve  B  the  pressure  taps  were 
still  in  the  flow  line  but  additional  pressure  taps  had  been 
drilled  along  the  length  of  the  core  and  curve  C  shows  the  result 
when  the  pressure  taps  were  moved  from  the  flow  line  to  the  end 
plates  of  the  mounted  core.  A  significant  change  in  the 
measured  flow  characteristics  is  apparent.  The  curves  for 
run  A  and  B  in  figure  8  indicate  that  the  greater  part  of  the  flow 
was  in  the  transition  region  while  curve  C  indicates  almost 
continuous  laminar  flow  for  the  range  of  flow  rates  covered. 

The  plot  of  the  pressure  squared  along  the  length  of  the  core, 
shown  in  figure  6 verifies  that  laminar  flow  existed  through 
the  range  of  flow  rates  covered. 

From  the  results  of  the  experimental  data  obtained  on 
the  long  core  samples  it  is  apparent  that  it  is  necessary  to 
measure  the  flowing  upstream  and  downstream  pressures  at  the 
core  face.  While  in  liquid  flow  tests  the  error  due  to  measuring 
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the  pressures  in  the  flow  lines  may  not  be  significant  when  the 
flow  rate  is  low,  when  the  flowing  fluid  is  gas  the  error  due  to 
the  inclusion  of  short  lengths  of  flow  line,  even  at  low  flow 
rates,  can  be  significant.  Figure  7  is  ample  proof  of  this  fact 
as  it  shows  the  considerable  change  in  the  experimental  flow 
characteristics  due  to  the  pressure  measuring  point. 

Flow  data  were  obtained  on  14  small  core  samples  using  the 
core  holder  apparatus  shown  in  figure  1.  Tables  5  through  36 
in  Appendix  I  summarize  the  measured  and  calculated  results  of 
the  experimental  runs.  The  calculated  flow  data  results  are 
shown  plotted  in  Figures  7  through  25.  The  values  of  o(  and  p  are 
obtained  from  the  flow  data  plotted  in  accordance  with  Equation 
5  and  show  the  straight  line  result  predicted  by  that  equation. 
Figure  9  shows  flow  data  for  the  three  cardium  sandstone  cores. 

At  low  mean  pressures  all  of  the  plotted  data  showed  deviations 
due  to  slippage  or  molecular  streaming.  The  low  permeability 
cores  showed  the  more  pronounced  deviations  due  to  the  phenomena 
as  shown  in  Figure  10  for  Sturgeon  Lake  limestone  cores,  and 
Figure  11  for  Viking  sandstone  cores. 

The  Indiana  limestone  is  a  low  permeability  rock  and 
exhibits  the  phenomena  of  gas  slippage  until  high  mean  pressures 
are  reached.  This  was  discussed  previously  in  relation  to  the 
compared  results  of  the  flow  data  obtained  on  the  long  and  short 
core  samples  of  the  Indiana  limestone.  When  low  permeability 
core  samples  are  being  used  to  obtain  flow  characteristics  for 
the  particular  rock, it  is  necessary  to  attain  rates  of  flow, 
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or  mean  pressures ,  where  the  effect  of  gas  slippage  is  not 

significant.  The  effect  of  gas  slippage  can  be  held  constant 

by  obtaining  data  at  a  constant  mean  pressure.  This  approach 

( 29) 

was  used  by  Stewart  and  Owens  in  a  study  of  turbulent 

and  laminar  flow  in  heterogeneous  limestones.  Flow  data  plots 
where  different  constant  mean  pressures  are  used  should  show 
a  series  of  curves  which  become  asymptotic  to  the  straight  line 
covering  the  flow  region  where  the  effect  of  slippage  is  negligible. 

Table  3  summarized  the  porosity  and  flow  data  obtained 
from  the  experimental  results. 


Table  3 

SUMMARY  OF  POROSITY  AND  FLOW  DATA  FOR  CONSOLIDATED  POROUS  MEDIA 


Core  Porosity  {%)  o <  (ft“8)  /3  (ft~~*~)  K  (md) 


1 

19.50 

4.30 

X 

„  3-2 

10  -i  r 
10lQ 

1.805 

X 

10° 

lo14 

21.89 

2 

2.36 

1.922 

X 

1.07 

X 

0.0049 

3 

18.82 

5.665 

X 

10 12 
lo12 

2.533 

X 

108 

108 

16.61 

4 

21.11 

3.18 

X 

1.85 

X 

29.59 

5 

17.16 

5.24 

X 

1072 

2.165 

X 

108 

17.96 

6 

14.55 

18.54 

X 

1013 

2.613 

X 

10 10 

0.508 

7 

13.60 

7.125 

X 

10*3 

1.807 

X 

1010 

1.32 

8 

19.36 

5.745 

X 

1011 

1.144 

X 

107 

163.81 

9 

19.46 

3.463 

X 

10^3 

1.972 

X 

107 

2.715 

10 

19.96 

3.85 

X 

ioiJ- 

1.107 

X 

107 

244.44 

11 

20.39 

2.265 

X 

1011 

5.689 

X 

106 

-1  *1 

415.50 

12 

9.82 

4.778 

X 

1014 

1.962 

X 

10  ^ 

0.197 

13 

10.28 

1.64 

X 

10^4 

3.417 

X 

10p 

0.574 

14 

18.85 

1.568 

X 

1014 

1.264 

X 

1010 

0.6002 
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Figures  3  and  12  show  the  plotted  flow  data  of  flow  rate 
against  the  difference  of  the  squared  upstream  and  downstream 
flowing  pressures.  The  higher  permeability  cores  gave  a  laminar 
flow  result  at  low  flow  rates  but  a  gradual  departure  from  the 
linear  relationship  was  soon  apparent.  This  gradual  transition 
zone  was  the  flow  regime  in  which  all  the  remainder  of  the  flow 
data  were  obtained.  The  expected  fully  developed  turbulent  region 
in  which  the  slope  of  the  plotted  results  would  have  a  constant 
value  of  2  was  not  achieved.  For  the  low  permeability  cores  the 
pressure  necessary  to  attain  that  region  would  have  been  far 
beyond  the  safety  limits  of  the  equipment,  while  for  the  high 
permeability  cores  the  flow  rate  necessary  would  be  greater 
than  the  rated  capacity  of  the  flow  meters  which  were  used. 
Certainly  the  plotted  curves  were  still  in  the  transition  region, 
however  a  greater  coverage  of  the  flow  regions  could  be  obtained 
with  a  large  gas  supply  and  high  flow  rate  measuring  equipment. 

A  common  method  of  plotting  flow  data  is  the  Reynolds 
number  -  friction  factor  correlation.  Figures  13  and  14  show 
the  flow  data  plotted  as  such  a  correlation,  with  the  Reynolds 

number  and  friction  factor  determined  from  equations  (8)  and  (9)  . 

(9) 

The  work  of  Cornell  and  Katz  contains  a  similar  plot  with  the 
results  of  Green  and  Duwez^4^  for  flow  through  porous  metals 
included.  Friction  factor  -  Reynolds  number  data  for  four  sand¬ 
stones  were  extracted  from  reference  (9)  and  plotted  on  figure 
13.  It  can  be  seen  that  all  the  data  are  consistent,  falling  on 
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FIGURE  13  REYNOLDS  NUMBER  -  FRICTION  FACTOR  PLOT 
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FIGURE  14  REYNOLDS  NUMBER  -  FRICTION  FACTOR  PLOT 


a  smooth  curve  defined  by  the  equation  f  =  +  2.  The  forms 

KQ 

of  the  friction  factor  and  Reynolds  number  used  for  the  correlation 
plots  presented  in  this  paper  are  useful  mainly  for  correlating 
data.  The  way  in  which  this  type  of  plot  can  be  used  in  the 
solution  of  flow  problems  is  explained  in  reference  30. 

Table  4  shows  a  comparison  of  the  permeability  determined 
from  both  the  Klinkenberg  and  flow  data  plots. 


Table  4 


COMPARISON  OF  PERMEABILITIES  DETERMINED 
FROM  BOTH  FLOW  DATA  PLOTS  AND  KLINKENBERG  PLOTS 

Permeability  (millidarcies) 


Core 


Flow  Data 
Plot 


Klinkenberg 

Plot 


% 

Deviation 


Core 

Core 

Core 

Core 

Core 


1  Cardium  Sandstone 
6  Sturgeon  Lake 
Limestone 

8  Berea  Sandstone 

9  Indiana  Limestone 

14  Belly  River  Sandstone 


21.89 

0.508 

163.8 

2.715 

0.600 


20.2 

0.557 

179 

2.68 

0 . 643 


-  7.72 

+  11.22 

+  9.28 

-  1.29 

-I-  7.17 


The  Klinkenberg  permeability  depends  on  the  flow  being  in  the 
laminar  region.  In  the  case  of  core  1,  the  Cardium  sandstone, 
figure  16  shows  that  the  flow  was,  at  best,  just  bordering  on 
the  laminar  flow  region.  If  it  were  still  in  the  transition  region, 
as  figure  8  appears  to  indicate,  the  permeability  determined 
would  be  too  low.  It  is  considered  likely  that  had  lower 
mean  pressures  been  reached  the  permeability  would  have  been 
higher.  The  flow  data  permeability  determination  depends  on  the 
elimination  of  the  effect  of  slippage,  an  effect  which  is  not 
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FIGURE  15  KLINKENBERG  PERMEABILITY  PLOT 


FIGURE  16  KLINKENBERG  PERMEABILITY  PLOT 
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accounted  for  in  flow  equation  (5).  If  what  is  assumed  to  be  the 
straight  line  region,  on  the  flow  data  plots  to  obtain  of  and 
P  ,  is  in  fact  still  slightly  inf luenced  by  slippage,  the  flow 
data  curve  would  in  fact  have  slight  curvature  up  to  higher  mean 
pressures.  The  visual  fit  of  a  straight  line  to  the  data  would 
give  an  intercept  Pf  which  would  be  too  low.  This  would  result 
in  a  permeability  value  which  would  be  too  high.  While  the  values 
of  permeability  determined  are  in  quite  good  agreement  the  fact 
that  the  flow  data  permeability  is  higher  in  most  cases  than  the 
Klinkenberg  permeability  lends  weight  to  the  argument  that  data 
runs  should  be  carried  out  for  at  least  two  constant  mean  pressures 
to  ensure  that  the  straight  line  region  of  the  plot  is  in  fact 
not  influenced  by  slippage. 

A  primary  aim  of  this  investigation  was  to  examine  the 
result  of  additional  flow  data  determined  characteristics  when 
correlated  on  a  plot  of  against  K  .  The  data  obtained  during 
this  investigation  together  with  that  of  Cornell  and  Katz^^  is 
shown  plotted  on  figure  17.  A  similar  plot  was  used  by  Cornell 
and  Katz  to  show  that  an  approximate  straight  line  relationship 
results.  The  equation  of  the  straight  line  was  found  to  be 


=  4.11  x  1010  (XI) 

k4/3 

(32)  e 

Subsequently  Janicek  and  Katz  recorrelated  the  data  of 

Cornell  and  Katz^"^  in  the  following  manner.  A  plot  of  the  formation 

factor  against  1/0  was  interpreted  to  give  the  linear  relationship 
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FIGURE  17  CORRELATION  OF  (3  FACTOR  WITH  PERMEABILITY  AND  POROSITY 
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F  =  4/0 


(12) 


A  comparison  of  the  equation  proposed  by  Green  and  Duwez 


(14) 


with  that  of  Cornell  and  Katz  led  to  the  relationship  that 


o( 


32  F 

ki  de: 


and  B  =  32  F  / 

'  *2  DE  ^T75 


which  with  the  substitution  of  k^  =  0.5,  F  =  4/0  and  K 


13 

9.411  x  10  /(X  results  in  the  equation 


D 


E  = 


8  (FK) 


1/2 


9.7  x  106 


(13) 


Observing  that  k2  is  related  to  porosity  and  permeability, 
k2  was  plotted  against  K/0  and  a  straight  line  relationship  was 


assumed  to  fit  the  data,  the  equation  of  which  was 

~)3/4 


k2  =  2.82  x  10 


-2 


(14) 


The  final  substitution  of  the  expressions  for  DE,  and  F  in 
the  relationship  expressed  above  led  to  the  equation  for 
in  terms  of  the  permeability  and  porosity 


5.5  x  10 

5/4  .3/4 

K  0 


(15) 


which  is  the  basis  of  the  constant  porosity  lines  of  Janicek  and 
Katz  plotted  on  figure  17. 

The  porosity  lines  so  developed  do  not  appear  to  fit  the 
experimental  data  points  with  any  consistency.  They  appear  to  be  too 
close  together.  For  this  reason  the  data  of  Cornell  and  Katz, 
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for  the  factors  k 2*  K,  0  and  F,  was  replotted  during  this  study 
as  k2  against  combinations  of  the  other  three  factors.  It  was 
hoped  that  a  grouping,  different  from  that  used  by  Janicek  and 
Katz,  would  lead  to  a  more  consistent  set  of  porosity  lines. 

One  of  the  groupings  used  was  a  dimensionless  form  where  k2  was 
plotted  against  0/F.  Again  a  straight  line  can  be  fitted  to 
the  data  points.  An  equation  for  k2  can  be  determined  and  an 
expression  for  ,  similar  to  that  above,  obtained.  This 
dimensionless  relationship  lends  considerable  weight  to  the 
porosity  and  very  little  to  the  permeability,  while  the  relation¬ 
ship  developed  by  Janicek  and  Katz  appears  to  weigh  the  effect 
of  the  permeability  very  strongly  and  minimizes  the  effect  of 
the  porosity.  Other  forms  of  plotting  the  data  examined  were 
k2  against  K0,  V  K0  and  K0/F.  These  plots  all  weigh  the  effect 
of  the  permeability  and  the  porosity  to  varying  degrees.  The 
form  of  the  plotted  data  which  was  used  was  that  of  k2  against 
K0/F  resulting  in  the  equation 


k 


2 


5.62 


(16) 


which  when  substituted  with  the  expression  for  and  F  results 
in  the  equation 

7 

_  6.6  x  10  ( 17 ) 

"  K1-13  0  2-76 

This  equation  is  the  basis  of  the  set  of  porosity  lines  on  figure 


17a. 
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FIGURE  17a  CORRELATION  OF  /3  FACTOR  WITH  PERMEABILITY  AND  POROSITY 


The  porosity  lines  so  determined  appear  to  fit  the 
experimental  data  with  much  more  consistency  than  those  predicted 
previously.  However  the  result  is  purely  empirical  and  can  only 
be  justified  if  shown  to  be  consistent  with  experimental  flow 
characteristic  data.  It  is  overly  optimistic  to  assume  that 
the  three  factors,  inertial  resistance  coefficient  p>  ,  permeability 
and  porosity  should  completely  characterize  the  flow  channel 
geometry  of  porous  media.  The  investigation  of  other  measurable 
properties  of  porous  media,  for  their  modifying  influence  on 
the  values  of  K,  0,  and  ,  is  necessary.  Some  of  those  which 
could  be  investigated  are,  the  fraction  of  dead  end  pore  volume, 
the  surface  roughness,  pore  size  distribution,  and  tortuosity. 

At  present  it  can  be  generally  said  that  the  factor  increases 
with  decreasing  permeability  and  decreasing  porosity* 

The  effect  of  the  temperature  on  the  porous  media  flow 

/  /,  (15) 

factors  c\  and  6  was  investigated  by  Greenberg  and  Weger 

for  the  flow  of  gas  through  porous  sintered  metals.  They 

found  that  as  the  temperature  was  increased,  the  flow  data 

curve  was  displaced  to  give  a  larger  value  of  ,  that  is,  the 

permeability  was  decreased.  The  ^  factor  however  appeared  to  be 

unchanged  as  the  slope  of  the  curves  remained  constant.  Johnson 

(19) 


and  Taliaferro 


conducted  controlled  temperature  flow  tests 


at  a  series  of  temperatures  with  the  stated  result  that  for  the 
test  procedure  used  no  apparent  change  in  the  permeability  with 
temperature  was  indicated.  Flow  data  at  a  series  of  temperatures 


•Xlv 


for  flow  through  a  Cardium  sandstone  sample  were  obtained  and 
are  shown  plotted  on  figure  18.  The  plotted  data  appear  to  verify 
the  results  of  Greenberg  and  Weger  with  they3  factor  unchanged 
and  a  decreasing  permeability  with  increased  temperature. 

The  factor  ^  is  a  measure  of  the  energy  loss  due  to  the 
inertial  effects  becoming  significant.  Due  to  the  complicated 
nature  of  the  flow  path,  the  effect  will  become  significant  at 
low  velocities  of  flow.  With  increased  temperature  very  little 
change  would  be  expected  in  the  flow  path  structure,  certainly 
not  enough  to  cause  an  increase  in  the  energy  loss  due  to  the 
inertial  effects.  It  is  reasonable  and  expected  that  in  the  case 
of  the  {2  factor  no  temperature  effect  would  be  apparent. 

The  permeability,  which  is  the  inverse  of  the  factor  of  , 
is  a  viscous  resistance  coefficient  and  as  such  is  dependent  on 
the  dimensions  of  the  flow  channels.  With  an  increased  temperature 
a  change  in  the  flow  channel  dimensions  might  be  expected. 

The  core  holder  used  was  a  steel  holder  with  the  core  mounted 
in  a  rubber  stopper.  The  stopper  was  held  in  the  tapered  section 
of  the  holder  as  shown  in  figure  1  with  the  follower  holding  the 
rubber  stopper  in  the  tapered  section  as  tightly  as  possible. 

The  change  in  the  flow  characteristic  <?(  of  the  core  should  not 
have  been  due  to  the  method  of  holding  the  core  but  must  have 
been  due  to  rock  expansion  alone.  The  result  of  the  increased 
temperature  would  be  an  expansion  of  the  rock  causing  a  reduction 
in  the  cross  sectional  area  of  the  flow  channel.  This  would 
cause  an  increase  in  the  flow  resistance  and  result  in  a  reduced 
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FIGURE  18  FLOW  DATA  PLOT  TO  OBTAIN  <*  AND  /3 


permeability,  which  is  what  was  experienced.  Greenberg  and 
Weger  state  that  since  the  expansion  of  the  metal  composing  the 
sintered  cores  is  negligible, the  reduction  in  permeability  must 
be  due  to  the  surface  irregularities,  particularly  where  the 
flow  passages  are  narrow.  They  feel  that  the  major  part  of  the 
viscous  resistance  takes  place  at  these  narrow  openings  and 
any  slight  change  in  the  size  or  position  of  the  irregularities 
v/ill  have  a  very  large  effect.  This  explanation  is  felt  to  be 
justified,  however  there  are  undoubtedly  many  such  narrow  points 
in  the  flow  paths  through  a  porous  rock  which  will  be  affected 
by  the  slight  rock  expansion  when  the  temperature  rises.  The 
slight  decrease  in  permeability  experienced  over  the  temperature 
range  covered  could  be  due  to  the  slight  decrease  in  cross 
sectional  area  at  these  constrictions / causing  the  increased 
viscous  flow  resistance. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

The  plot  of  yd  against  K  with  the  parameter  of  porosity 
shows  that  there  appears  to  be  a  relationship  relating  these 
variables.  The  manner  in  which  this  relationship  may  be 
completely  expressed  is  not  apparent  other  than  by  the  empirical 
relationship  which  has  been  developed.  The  effects  of  other 
factors  on  the  three  variables  need  to  be  examined  to  try  to 
find  a  possible  relationship  between  the  variables  which  is 
applicable  to  all  porous  media.  Until  such  a  relationship  is 
developed  the  modified  expression  for  ,  equation  (17) , 
together  with  the  friction  factor  -  Reynolds  number  correlation 
can  be  used  for  gas  flow  calculations  when  the  permeability 
and  porosity  are  known. 

The  effect  of  gas  slippage  is  considered  to  be  of 
significance  only  at  low  mean  pressures;  however  the  flow 
data  plots  to  obtain  the  inertial  resistance  coefficient  j9  and 
permeability  K  ,  depend  on  the  elimination  of  this  effect  to 
predict  the  correct  values  of  the  two  factors.  Several  flow 
tests  for  each  core,  in  which  experimental  data  were  obtained 
at  different  constant  mean  pressures,  would  show  the  flow  rate 
up  to  which  the  slippage  effect  was  significant.  It  is  re¬ 
commended  that  flow  data  on  cores,  covering  a  wide  range  of 
permeability,  be  obtained  at  different  constant  mean  pressures. 
The  plots  of  these  data  would  indicate  the  flow  rates  that 
probably  have  to  be  achieved,  for  any  permeability  range,  before 
the  assumption  that  the  data  points  fit  the  straight  line, 
predicted  by  the  flow  equation  (5),  is  necessarily  true. 
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A  further  investigation  on  the  effect  of  temperature  on 
the  flow  characteristics  of  porous  rocks  needs  to  be  carried 
out.  Changes  in  the  experimental  apparatus  setup  would  allow 
better  control  of  the  temperature.  A  core  holder  which  could 
be  immersed  in  a  constant  temperature  bath  should  enable 
positive  control  of  the  temperature  to  be  maintained.  With 
a  core  holder  such  as  the  one  used  in  this  study  it  is 
necessary  that  the  follower  be  kept  tightly  screwed  down. 

This  has  to  be  checked  on  every  change  of  temperature. 

The  small  core  samples  used  in  this  study  were  taken 
from  bulk  samples  of  the  formation  rock.  In  many  cases  small 
cores  samples  were  obtained  from  immediately  adjoining  sections 
of  the  bulk  sample,  while  the  flow  characteristic  factors 
determined  for  each  small  core  sample  are  significant,  they 
have  to  be  interpreted  as  applying  to  the  core  on  which  they 
were  determined,  thus  different  permeability  and  ft  factor  values 
were  found  for  each  small  cardium  core  studied.  These  factors 
describe  the  properties  of  the  particular  core,  and  no  one  small 
core  can  be  said  to  characterize  cardium  sandstone  generally. 
Characteristic  factors  based  on  the  average  properties  of  a 
cardium  sandstone  reservoir  would  however  be  significant  as 
characterizing  the  bulk  rock. 

The  difficulty  that  was  experienced  in  mounting  large 
core  samples  has  led  to  a  continuing  investigation  into 
different  mounting  techniques.  The  use  of  an  epoxy  resin 
material  for  mounting  large  core  samples  has  been  modified,  to 
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eliminate  the  penetration  of  the  core  by  the  epoxy  resin,  which 
is  highly  wetting.  While  the  use  of  this  mounting  material  is 
of  doubtful  value  for  large  core  samples  it  can  be  used  success^ 
fully  for  small  cores.  The  methods  which  are  felt  to  hold  the 
most  promise  are  those  in  which  the  core  is  thinly  coated  with 
a  non  wetting  material  and  then  mounted  using  a  non  shrink 
material  or  other  sufficiently  strong  matieral. 

The  long  core  data  also  points  out  the  need  for  properly 
designed  equipment,  with  the  proper  location  of  pressure  taps 
and  other  measuring  points.  The  method  of  drilling  pressure  taps 
along  the  length  of  the  core  needs  investigation.  The  possibility 
of  pressure  taps  being  positioned  on  the  core  prior  to  its 
mounting  might  be  one  way  to  avoid  drilling  through  the  mounting 
material  and  into  the  core. 

The  Reynolds  number  -  friction  factor  correlation  shown 
in  figures  13  and  14  is  based  on  the  definitions  expressed  by 
equations  (8)  and  (9) .  With  the  Reynolds  number  so  defined,  the 
correlation  graphs  indicate  that  for  the  flow  of  a  gas  through 
porous  media,  the  deviation  from  purely  viscous  flow  takes  place 
at  the  same  Reynolds  number.  This  Reynolds  number  value  is  between 
0.08  and  0.10.  If  permeability  and  porosity  data  are  available 
it  is  possible  to  determine  the  mass  flow  rate  at  which  Darcy’s 
law  for  gas  flow  can  no  longer  be  used  to  obtain  calculated  flow 
data.  The  data  below  a  Reynolds  number  value  of  approximately  0.03 
appear  to  be  influenced  by  gas  slippage  to  a  varying  degree.  This 
is  indicated  by  the  data  points  falling  below  the  line  defined  by 
the  equation  f  =  2/NRe. 
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NOMENCLATURE 

coefficient  of  inertial  resistance  term  of  flow 
equation  (feet-1) 

difference 

pressure  (pounds  force/square  inch  absolute) 
mean  pressure  (pounds  force/square  inch  absolute) 

-  molecular  weight  of  a  gas 

conversion  factor  (32.17  pounds  mass  x  foot/pounds 
force  x  seconds  squared) 

length  of  porous  sample  (inches) 

compressibility  factor  of  gas 

-  gas  constant  (1543  pounds  force  x  cubic  feet/square 
feet  x  pound  mole  x  °R) 

-  absolute  temperature  °R 

viscosity  of  gas  (pounds  mass/foot  second) 
mass  velocity  (pounds  mass/square  foot  second) 

-  coefficient  of  laminar  term  of  flow  equation  (feet~^) 
permeability  (millidarcies) 

apparent  permeability  (millidarcies) 
porosity 

velocity  (feet/second) 
flow  rate  (cubic  feet/hour) 

cross  sectional  area  of  porous  sample  (square  inches) 
density  (pounds  mass/cubic  feet) 

Electrical  resistivity  factor;  electrical  resistance 
of  sample  saturated  with  brine/resistance  of  quantity 
of  brine  of  same  size  and  shape  as  bulk  volume  of  sample. 
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1 


k 
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D 


E 


N. 


Re 


dimensionless  geometrical  constant  in  laminar  term 
of  flow  equation 

dimensionless  geometrical  constant  in  inertial  term 
of  flow  equation 

effective  diameter  of  pore  structure  (feet) 

Reynolds  number  (dimensionless) 
friction  factor 
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TABLE  6(A) 

GAS  FLOW  DATA 

ALUNDUM  CORE  (TA  7921) 
(before  pressure  taps  drilled) 


Run 

No 

Pressure-Psia 

P1  p2 

Tf  low 
(average) 
°R 

Qc 

(SCFH) 

(P12-P22) 

(Psia)2 

G 

(lbs/f t^sec) 

1 

pavg 

(Psi”'*-) 

Ka 

(md . ) 

1 

40.4 

17.3 

534.7 

6.85 

1340 

.006452 

.0347 

650.79 

2 

54.9 

22.8 

535.7 

11.65 

2500 

.01097 

.0257 

594.18 

3 

65.8 

28.2 

537.4 

15.45 

3540 

.01455 

.0213 

558.91 

4 

75.4 

33.3 

538.9 

18.73 

4580 

.01764 

.0184 

524.61 

5 

82.7 

37.8 

538.9 

21.41 

5410 

.02017 

.0166 

508.08 

69 


TABLE  6(B) 

GAS  FLOW  DATA 

ALUNDUM  CORE  (TA  7921) 

(Pressure  Taps  Drilled  Along  Length  of  Core. 
Pressure  Taps  Upstream  and  Downstream  of  End  Plates.) 


Pressure 

-Psia 

^■f  low 

Qc 

(P].2-P22) 

G 

1 

Run 

(average) 

Pavg, 

Ka 

No 

Pi 

P2 

°R 

(SCFH) 

(Psia)z 

(lbs/ ft 

^sec) 

(Psi-1) 

(Md.) 

1 

213.1 

138.8 

531.5 

84.60 

26130 

7.969  x 

CM 

1 

O 

r— 1 $ 

0.00568 

419.8 

2 

212.5 

156.9 

530.8 

71.58 

20540 

6.742 

h 

0.00541 

442.43 

3 

213.7 

189.5 

530.8 

40.75 

9760 

3.838 

ft 

0.00496 

530.86 

4 

212.5 

173.6 

530,9 

56.90 

15040 

5.360 

11 

0.00518 

480.72 

5 

214.2 

204.1 

530.9 

20.12 

4220 

1.895 

II 

0.00478 

605.94 

6 

214.1 

208.2 

530.9 

13.47 

2490 

1.175 

ii 

0.00474 

636.59 

7 

215.1 

211.8 

530.9 

6.53 

1380 

0.615 

If 

0.00468 

604.41 

8 

213.6 

212,0 

530.3 

3.24 

660 

0.305 

ft 

0.00470 

624.88 

9 

21.58 

13.85 

531.3 

1.722 

274 

0.162 

ii 

0.0564 

790.81 

10 

21.75 

14.19 

532.0 

1.700 

272 

0,160 

ii 

0.0557 

787.80 

11 

21.47 

14.54 

532.0 

1.560 

250 

0.147 

f  1 

0.0556 

787.71 

12 

21.60 

17.02 

532.7 

1.097 

177 

0.103 

it 

0.0518 

783.27 

13 

21.66 

21.34 

532.7 

0.0541 

14.10 

0.00500 

It 

0.0465 

702.90 

14 

29.39 

22.11 

531.7 

3.89 

643 

0.366 

ii 

0.0452 

761.2 

15 

29.41 

26.11 

531.7 

2.18 

344 

0.205 

it 

0.0383 

799.43 

16 

28.49 

27.63 

532.0 

0.582 

94.5 

0.0548 

II 

0.0362 

776.09 

17 

37.26 

16.75 

529.2 

6.48 

1108 

0.610 

ii 

0.037 

731.47 

18 

36.51 

21.84 

529.2 

5.09 

856 

0.479 

ii 

0.0343 

743.22 

19 

36.22 

30.09 

529.2 

2.47 

407 

0.233 

II 

0.0302 

759.30 

20 

79.1 

38.89 

530.2 

22.73 

4570 

2.141 

1? 

0.0169 

600.77 

21 

79.4 

56.7 

530.2 

15.78 

3095 

1.486 

It 

0.0158 

641.38 

22 

78.9 

67.4 

530.9 

8.82 

1680 

0.831 

II 

0.0137 

662.76 

23 

79.3 

77.3 

530.9 

2.15 

313 

0.203 

ti 

0.0128 

865.93 

24 

113.3 

107.6 

530.9 

6.49 

1250 

0.611 

If 

0.00905 

656.3 

25 

114.2 

90.3 

530.9 

15.98 

3190 

1.505 

II 

0.00937 

632.71 

26 

114.5 

72.9 

530.9 

34.22 

7800 

3.223 

II 

0.0107 

553.44 

27 

114.3 

63.9 

530.9 

38.31 

8980 

3.609 

II 

0.0122 

538.39 

28 

163.8 

101.2 

530.9 

60.92 

16590 

5.738 

II 

0.00755 

463.93 

29 

164.2 

118.0 

529.7 

51.15 

13030 

4.818 

II 

0.0071 

495.67 

30 

163.7 

136.8 

529.7 

34.96 

8080 

3.293 

it 

0.0067 

547.03 

31 

163.3 

156.8 

530.2 

10.29 

2060 

0.969 

II 

0.0062 

630.81 

■ 

•  2 

Run 

No 

1 

2 

3 

4 

5 

6 

7 

3 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 
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TABLE  6 (Cl) 
GAS  FLOW  DATA 


ALUMDUN  CORE  (TA  7921) 
(Pressure  Taps  in  End  Plates) 


Pressure- 

-Psia 

Tf  low 

Qc 

(Pi2-p,2) 

G 

1 

(Average) 

o 

Pavg 

(Psi'1) 

Ka 

pi 

P2 

°R 

(SCF) 

(Psia) ^ 

(lbs/ ftzsec) 

(Hd.) 

63.5 

33.1 

535 

16.66 

2940 

1.569  x  10“2 

.0207 

722.44 

79.0 

45.5 

535 

23.24 

4170 

2.189 

.0161 

711.37 

94.34 

57.7 

535 

31.66 

5570 

2.982 

.0142 

726.53 

79.2 

45.6 

535 

24.05 

4180 

2.265  " 

.0160 

734.34 

20.40 

13.77 

536.3 

1.392 

227 

0.131  " 

.0585 

784.19 

23.21 

14.07 

536.3 

2.083 

341 

0.196  " 

.0536 

780.81 

25.44 

14.39 

536.3 

2.68 

440 

0.247 

.0502 

778.12 

18.02 

13.60 

536.5 

0.86 

140 

0.081 

.0633 

785.67 

33.03 

13.17 

536.5 

4.98 

818 

0.469 

.0433 

778.71 

40.80 

19.24 

536.5 

7.70 

1295 

0.725  " 

.0333 

760.38 

51.6 

24.81 

536.5 

11.61 

2045 

1.094 

.0262 

726.23 

14.43 

13.29 

535.4 

0.194 

31.5 

0.0183 

.0722 

784.57 

13.90 

13.31 

535.4 

0.10 

16.1 

0.0094 

.0735 

792.82 

13.48 

13.31 

535.4 

0.0326 

4.53 

0.0031  " 

.0747 

919.06 

77.8 

44.4 

536.3 

23.37 

4080 

2.201 

.0163 

733.66 

94.1 

58.5 

536.3 

30.92 

5440 

2.913  " 

.0131 

728.78 

115.9 

78.6 

536.1 

40.42 

7250 

3.807 

.0115 

714.27 

140.3 

101.4 

535.3 

53.16 

9400 

5.007 

.0083 

725.73 

169.8 

130.00 

534.4 

67.94 

11940 

6.400 

.0067 

730.79 

190.2 

148.3 

534.4 

80.23 

14170 

7.557 

.0059 

728.00 

214.4 

171.5 

534.1 

93.48 

16570 

8.805  " 

.0052 

726.37 

238.9 

194.7 

533.5 

106.2 

19190 

10.00 

.0046 

712.85 

256.44 

212.3 

533.5 

116.2 

20680 

10.94 

.0043 

723.42 

226.0 

171.1 

533.5 

121.6 

21795 

11.45  " 

0.005 

717.67 

195.3 

142.5 

533.5 

100.9 

17820 

9.501 

0.006 

726.48 

252.4 

196.1 

533.5 

139.7 

25270 

13.16 

0.0045 

712.42 

146.2 

97.6 

533.5 

68.25 

11860 

6.429 

0.008 

735.62 

Run 

No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 


-  71 


TABLE  6 (C 2 ) 

GAS  FLOW  DATA 

ALUNDUM  CORE  (TA  7921) 
(pressure  taps  in  end  plates) 


Pressure 

pl 

.-Psia 

P2 

If  low 
(average) 
°R 

Sc 

(SCFH) 

(Pl2-P22> 

(Psia)3 

G 

(lbs  ft 

2 

sec) 

1 

Pavg 

(psi"1) 

Ka 

(md . ) 

27.11 

13.49 

530 

3.466 

553 

3.264 

x  10"3 

0.0493 

784.0 

33.34 

14.06 

530 

5.719 

914 

5.386 

fl 

0.0422 

783.2 

39.3 

14.7 

530 

8.230 

1329 

7.753 

fl 

0.0370 

775.4 

59.6 

19.83 

529 

18.566 

3164 

17.49 

II 

0.0251 

729.1 

83.6 

32.50 

529 

34.99 

5933 

32.97 

it 

0.0172 

739.2 

97.1 

41.5 

529 

45.22 

7704 

42.58 

If 

0.0144 

724.5 

114.5 

54.0 

529 

59.07 

10180 

55.63 

IV 

0.0119 

727.9 

133.3 

67.9 

529 

74.57 

13160 

70.25 

tl 

0.0099 

712.2 

153.3 

83.6 

529 

92.48 

16500 

87.12 

f  1 

0.0084 

705.6 

173.7 

99.5 

528 

110.3 

20280 

103.9 

II 

0.0073 

684.2 

164.2 

92.6 

528 

100.3 

18380 

94.49 

II 

0.0078 

686.1 

143.8 

75.0 

529 

81.78 

15050 

77.03 

II 

0.0091 

683.2 

(£0)d  'UCAT 


- 

Run 

No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 
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TABLE  7 
GAS  FLOW  DATA 
CORE  #1  CARDIUM  SANDSTONE 


Pressure 

-Psia 

low 

Qc 

(Pi2-P22) 

G 

1 

(average) 

Pavg 

Pi 

P2 

°R 

(SCFH) 

(Psiar 

(Lbs  /£ t^sec) 

(Psi"1) 

66 . 1 

13.96 

534.2 

2.12 

4173 

0.804  x  10"2 

0.025 

80.1 

14.40 

534.4 

3.08 

6215 

1.167 

0.0212 

93.8 

15.04 

534.6 

4.18 

8571 

1.584 

0.0184 

115.9 

16.60 

534.6 

6.17 

13160 

2.339 

0.0151 

138.1 

18.95 

534.4 

8.48 

18720 

3.214 

0.0127 

153.4 

20.97 

534.4 

10.18 

23100 

3.859 

0.0115 

174.2 

24.05 

534.4 

12.56 

29770 

4.761 

0.0101 

189.3 

26.65 

534.4 

14.38 

35120 

5.451 

0.0093 

204.7 

29.40 

534.4 

16.30 

41030 

6.178 

0.0085 

218.3 

32.00 

534.4 

17.90 

46640 

6.785 

0.0080 

233.8 

35.13 

534.4 

19.83 

53420 

7.516 

0.0074 

242.6 

36.98 

534.4 

20.96 

57510 

7.945  " 

0.0072 

234.6 

35.32 

530.7 

20.90 

53810 

7.922  " 

0.0074 

271.9 

43.3 

530.1 

26.17 

72050 

9.920 

0 . 0064 

310.0 

52.3 

530.4 

31.84 

93360 

12.07  " 

0.0055 

361.2 

64.5 

530.8 

39.59 

126330 

15.01 

0.0047 

412.9 

77.7 

530.8 

47.89 

164440 

18.15 

0.0041 

460.3 

90.0 

529.6 

55.55 

203770 

21.06 

0.0036 

513.5 

103.9 

529.2 

64.57 

252920 

24.78 

0.0032 

810.8 

185.6 

529.0 

116.3 

623000 

44.07 

0.0020 

754.0 

170.2 

527.8 

106.1 

539540 

40.22  " 

0.0022 

700.6 

155.5 

527.4 

97.20 

466650 

36.84 

0.0023 

641.9 

139.3 

527.3 

86.68 

392610 

32.86 

0.0026 

589.3 

122.2 

527.2 

76.00 

332300 

28.81 

0.00281 

527.3 

108.4 

527.2 

67.58 

266300 

25.62 

0.00315 

Ka 

(Md . ) 


23.03 

22.49 

22.15 

21.32 

20.60 

20.06 

19.22 

18.67 

18.13 

17.53 
16.95 
16 . 66 
17.55 
16.42 

15.44 
14.24 
13.26 
12.41 
11.65 

8.63 

9.05 

9.55 

10.10 

10.44 

11.54 


Run 

No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 


i  j 


TABLE  8 
GAS  FLOW  DATA 

CORE  #2  CARDIUM  SANDSTONE 


Pressure-Psia  Tflow  Qc  (P].2-P22)  G 


Pl 

P2 

R 

(SCFH) 

(Psia)^ 

(lbs/ft^ 

sec) 

103.6 

13.62 

532.8 

0.00312 

10550 

0.1181  x 

10'4 

133.4 

13.62 

532.8 

0.00348 

17200 

0.1316 

It 

163.6 

13.62 

532.8 

0.00488 

2656 

0.1846 

It 

234.0 

13.42 

529.2 

0.00768 

5455 

0.2912 

tt 

317.1 

13.43 

529.2 

0.0128 

100390 

0.4835 

If 

315.7 

13.01 

530.7 

0.0121 

99520 

0.4574 

It 

415.3 

13.01 

530.7 

0.0193 

172290 

0.7305 

tl 

521.8 

12.96 

530.7 

0.0287 

272130 

1.086 

ft 

618.9 

12.96 

530.7 

0.0387 

382890 

1.464 

tl 

745.5 

13.27 

530.2 

0.0541 

555560 

2.047 

II 

896.5 

13.29 

530.2 

0.0748 

803480 

2.828 

II 

1027.1 

13.30 

530.2 

0.0984 

1054680 

3.588 

II 

1151.8 

13.31 

530.8 

0.1153 

1326540 

4.370 

II 

1320.0 

13.33 

530.9 

0.1445 

1742090 

5 .466 

II 

1477.3 

13.35 

530.9 

0.1729 

2182210 

7.541 

II 

1575.0 

13.39 

530.9 

0.1964 

2480290 

7.430 

tt 

1723.4 

13.48 

531.2 

0.2203 

2969890 

8.334 

II 

1885.9 

13.48 

531.2 

0.2540 

3556440 

9.609 

II 

Run 

No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 
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TABLE  9 
GAS  FLOW  DATA 
CORE  #3  CARD I UN  SANDSTONE 


Pressure 

pi 

-Psia 

P2 

Tf  low 
(Average) 
°R 

Qc 

(SCFH) 

(Pi2-p22) 

(Psia)2 

G 

(lbs  /ft2sec) 

20.59 

13.54 

540.2 

0.100 

240.6 

.0379  x  10 

32.11 

13.55 

539.7 

0.340 

847.5 

.1286 

53.5 

13.60 

537.9 

1.024 

2672 

.3874 

18.1 

13.65 

537.9 

2.417 

6395 

.9144 

114.2 

13.79 

537.9 

4.594 

13850 

1.738 

163.3 

14.24 

537.9 

8.73 

26480 

3.303 

188.1 

14.59 

537.9 

11.02 

35150 

4.169 

239.3 

15.71 

537.9 

16.24 

57020 

6.144 

283.3 

17.15 

535.0 

20.93 

79960 

7.918 

338.2 

19.41 

535.0 

28.00 

114000 

10.59 

413.6 

23.80 

534.4 

37.41 

170500 

14 . 15 

488.5 

29.08 

533.1 

47.38 

237800 

17.93 

577.1 

36.03 

533.4 

59.32 

331800 

22.44 

642.4 

41.46 

533.2 

68.63 

410900 

25.96 

1216.7 

97.7 

527.9 

158.2 

1471000 

59.87 

1107.6 

85.8 

527.4 

139.4 

1219500 

52.72 

984.9 

74.3 

527.4 

120.4 

964600 

45.55 

843.7 

61.2 

527.4 

97.81 

708000 

37.00 

743.0 

52.1 

527.5 

83.46 

549300 

31.58 

CAT 


' 

.  . 

• 

Run 

No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 
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TABLE  10(A) 

GAS  FLOW  DATA 
CORE  #4  CARDIUM  SANDSTONE 
(Tf  =  73°F) 


Pressure-Psia 

PL  p2 

1  flow 
(average) 
°R 

Qc 

(SCFH) 

<P1  -P2  > 
(Psia)2 

G 

(lbs/ ft^sec) 

801.8 

88.4 

531.8 

142.1 

635100 

53.87  x  10"2 

736.1 

79.7 

530.9 

128.9 

535500 

48.87 

665.8 

70.3 

530.8 

113.8 

438400 

43.14 

601.3 

62.0 

530.85 

99.30 

357700 

37.64 

499.4 

48.6 

531.05 

77.93 

247000 

29.54  " 

430.6 

40.1 

531.85 

63 . 66 

183800 

24.13 

370.6 

32.76 

532.6 

51.84 

136300 

19.65  " 

314.7 

26.48 

533.35 

41.05 

98300 

15.56 

260.5 

21.21 

533.85 

30.94 

67400 

11.73 

206.3 

17.36 

537.9 

21.07 

42240 

7.987 

147.6 

14.84 

538.0 

12.47 

21560 

4.727 

101.9 

13.87 

538.0 

6.48 

10200 

2.456  " 

74.3 

13.61 

538.0 

3.54 

5335 

1.340 

44.4 

13.54 

538.0 

1.62 

1785 

0.614 

1 

Run 

No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 
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TABLE  10(B) 

GAS  FLOW  DATA 

CORE  ii 4  CARDIUM  SANDSTONE 
(Tf  *  95°F) 


Pressure 

pi 

-Psia 

P2 

^  f  1  ow 
(average) 
°R 

Qc 

(SCFH) 

(pl2-p22) 

( Psia) ^ 

G 

(lbc/ft2 

sec) 

86.1 

13.63 

554.25 

4.57 

7224 

1.732  x 

icf2 

56.4 

13.46 

552.2 

1.99 

3000 

0.7543 

II 

138.8 

14.42 

555.3 

10.74 

19070 

4.071 

It 

211.5 

17.24 

554.85 

21.08 

44400 

7.990 

fl 

301.4 

24.23 

554.95 

36.59 

90300 

13.87 

It 

364.2 

30.97 

554.2 

48 . 05 

131700 

18.21 

II 

415.9 

36.56 

555.0 

57.51 

171600 

21.80 

II 

458.9 

41.9 

554 . 6 

66.07 

208800 

25.04 

II 

510.0 

48.3 

555.65 

76.09 

257700 

28.84 

M 

603.8 

60.1 

555.65 

94.93 

360900 

35.98 

it 

809.6 

86.5 

555.1 

138.5 

647900 

52.49 

II 

721.1 

75.2 

556.3 

119.6 

514400 

45.35 

it 

656.9 

68.1 

555.05 

106.2 

426900 

40.27 

II 

-  77 


TABLE  10(C) 

GAS  FLOW  DATA 

CORE  #4  CARD I UK  SANDSTONE 


(Tf 


Run 

No 

Pressure-Psia 

?!  P2 

P  f  1  ow 
(average) 
°R 

1 

260.1 

20.23 

575.0 

2 

743.6 

76.6 

574.8 

3 

658.8 

66.0 

574.3 

4 

580.0 

55.90 

575.45 

5 

498.4 

45.9 

575.3 

6 

420.3 

36.6 

575.45 

7 

328.9 

26.3 

575.0 

8 

236.2 

18.5 

574.65 

=  115°F) 


(SCFH) 

(Pi2-f22) 

2 

(Psia) 

G 

2 

(lbs/ft  sec) 

27.76 

67230 

10.52  x  10‘2 

120.2 

546600 

45.57  " 

102.5 

429700 

38.84 

86.55 

333200 

32.81 

70.68 

246300 

26.79 

56.01 

175300 

21.23  " 

39.49 

107500 

14.97 

24.11 

55430 

9.14 

Run 

No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 
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TABLE  11 
GAS  FLOW  DATA 

CORE  #5  CARDIUM  SANDSTONE 


Pressure 

pi 

-Psia 

P2 

1  flow 
(average) 
°R 

<5c 

(SCFH) 

(pl  -p2  > 

( Psia)2 

G 

(lbs /ft 

2sec) 

1014.5 

86.5 

529.2 

140.5 

1021700 

53.36  x 

10"  2 

905.5 

75.4 

528.7 

122.30 

814300 

46.45 

n 

1096.1 

95.0 

525.4 

155.01 

1192400 

58.88 

it 

1087.1 

94.0 

525.4 

154.2 

1173000 

58.57 

ii 

809.8 

65.7 

526.8 

106.1 

651400 

40.30 

M 

694.0 

54.3 

530.3 

87.30 

478700 

33.16 

it 

616.2 

47.40 

530.1 

74.08 

378000 

28.14 

ti 

537.6 

39.21 

530.4 

62.24 

287400 

23.64 

II 

461.8 

32 . 14 

530.6 

50.67 

212200 

19.25 

II 

359.8 

23.53 

533.4 

35.64 

128900 

13.54 

II 

306.0 

19.96 

533.2 

28.16 

93250 

10.70 

II 

133.0 

14.04 

534.6 

7.56 

17500 

2.87 

IV 

213.1 

15.77 

534.7 

16.45 

45200 

6.25 

II 

97.3 

13.73 

535.7 

3.78 

9276 

1.44 

If 

60.5 

13.60 

536.7 

1.68 

3477 

0/64 

It 

25.16 

13.56 

533.8 

0.256 

449 

0.097 

II 

(a;  r  vr?) 

o.  m 

.0 

Run 

No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 
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TABLE  12 
GAS  FLOW  DATA 

CORE  #6  STURGEON  LAKE  SOUTH  LIMESTONE 


Pressure 

-Psia 

low 

Qc 

(P12-P22) 

G 

1 

(average) 

r\ 

Pavg 

Ka 

pi 

P2 

°R 

(SCFH) 

(Psia)z 

(lbs/ftz 

sec) 

(Psi"1 

(Md.) 

17.93 

13.61 

531.2 

0.00368 

136 

0.141  x 

10'4 

0.0634 

1.29 

22.91 

13.61 

531.7 

0.00849 

340 

0.324 

91 

0.0548 

1.19 

30.67 

13.61 

531.8 

0.0172 

756 

0.657 

ii 

0.0452 

1.08 

38.15 

13.61 

531.8 

0.0269 

1270 

1.028 

ii 

0.0386 

1.005 

48.00 

13.63 

528.9 

0.0411 

2120 

1.571 

19 

0.0324 

0.917 

55.10 

13.59 

530.2 

0.0536 

2846 

2.048 

it 

0.0291 

0.894 

64.00 

13.57 

530.7 

0.0711 

3914 

2.717 

ii 

0.0258 

0.864 

76.5 

13.57 

530.7 

0.0987 

5666 

3.771 

ii 

0.0221 

0.829 

89.7 

13.58 

530.7 

0.131 

7963 

5.002 

99 

0.0194 

0.783 

99.2 

13.59 

530.7 

0.158 

9650 

6.026 

ii 

0.0177 

0.780 

299.1 

13.72 

530.7 

1.084 

89260 

41.40 

99 

0.0064 

0.579 

398.3 

13.93 

530.7 

1.806 

158440 

69.00 

it 

0.0049 

0.549 

460.3 

14.09 

530.7 

2.221 

211720 

84.88 

ii 

0.0042 

0.507 

524.2 

14.29 

530.7 

2.749 

274570 

105.0 

ii 

0.0037 

0.485 

573.1 

14.49 

530.9 

3.189 

328190 

121.9 

ii 

0.0034 

0.472 

601.8 

14.61 

532.5 

3.439 

361890 

131.4 

ii 

0.0032 

0.462 

666.8 

14.91 

532.5 

4.06 

444450 

155.1 

ii 

0.0029 

mm 

711.3 

15.15 

532.5 

4.50 

505700 

171.95 

it 

0.0028 

- 

775.8 

15.59 

532.5 

5.17 

601600 

197.6 

ii 

0.0025 

- 

1456.8 

24.15 

532.5 

13.59 

2121800 

519.3 

it 

0.0014 

- 

1370.4 

22.68 

531.8 

12.39 

1877400 

473.4 

ii 

0.0015 

- 

1283.8 

21.30 

531.8 

11.23 

1647700 

429.10 

n 

0.0016 

- 

1198.3 

20.07 

530.2 

10.10 

1435500 

385.9 

ii 

0.0017 

- 

1123.0 

19.12 

530.9 

9.22 

1260700 

3.523 

ti 

0.00192 

- 

1006.6 

17.70 

530.9 

7.74 

1012900 

2.957 

99 

0.00195 

mm 

917.4 

16.76 

530.2 

6.68 

841300 

2.552 

ii 

0.00214 

- 

830.3 

15.98 

530.2 

5.68 

689100 

2.170 

91 

0.00236 

- 

. 

Run 

No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 


Pressure 

pl 

22.57 

39.98 

61.4 

92.7 

207.5 

308.9 

469,4 

604.3 
761.0 

910.4 
1060.7 

1202.4 
1473.1 
1390.6 
1306.0 
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TABLE  13 
GAS  FLOW  DATA 

CORE  #7  SOUTH  STURGEON  LAKE  LIMESTONE 


Psia 

P2 

Tf  low 
(average) 
°R 

Qc 

(SCFH) 

(pi2-p22) 

(Psia)2 

G 

(lbs/ft2 

sec) 

13.52 

537.2 

0.0187 

326.7 

0.7102  x 

10"4 

13.52 

537.3 

0.0744 

1416. 

2.823 

M 

13.52 

537.5 

0.172 

3583 

6 . 54 

t! 

13.52 

537.5 

0.374 

8413 

14.20 

ii 

13.57 

538.7 

1.487 

42850 

56.43 

tr 

13.62 

538.7 

2.797 

95250 

106.12 

f  i 

13.81 

538.7 

5.37 

220100 

203.76 

ii 

14.08 

538.7 

7.76 

365000 

294.44 

ii 

14.55 

538.2 

10.73 

578800 

407.13 

ft 

15.19 

537.2 

13.61 

828500 

516.41 

ft 

16.01 

536.3 

16.87 

1124800 

640.11 

it 

16.95 

535.95 

19.86 

1445500 

753.56 

ft 

19.35 

534.1 

26.83 

2169500 

1018.0 

it 

18.57 

533.2 

24.75 

1933300 

939.1 

tt 

17.84 

533.2 

22.28 

1705200 

845.4 

it 

. 

■ 

Run 

No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 
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TABLE  14 
GAS  FLOW  DATA 
CORE  #8  BEREA  SANDSTONE 


Pressure-Psia 

P1  p2 

Tf  low 
(average) 
°R 

Qc 

(SCFH) 

(Pt2-p22) 

(Psia)^ 

G 

(lbs/f t^sec) 

1 

Pavg 

(Psi’^) 

14.68 

13.75 

530.8 

0.133 

26.6 

0.051  x 

<N 

1 

O 

r— 4 

0.0704 

14.19 

13.75 

530.8 

0.0638 

12.43 

0.024 

II 

0.0716 

19.01 

13.85 

530.8 

0.834 

170 

0.317 

VI 

0.0609 

30.02 

14.56 

530.8 

3.26 

690 

1.239 

ti 

0.0449 

40.17 

16.19 

530.8 

6.11 

1352 

2.323 

VI 

0.0355 

63  •  6 

23.98 

526.3 

13.66 

3465 

5.192 

II 

0.0229 

83.3 

33.76 

527.0 

20.96 

5801 

7.965 

VI 

0.0171 

83.4 

33.68 

527.0 

20.82 

5818 

7.912 

IV 

0.0171 

114.7 

51.4 

527.0 

33.30 

10510 

12.65 

II 

0.0120 

133.2 

63.0 

529.1 

40.89 

13760 

15.54 

II 

0.0102 

153.5 

75.5 

529.1 

48.93 

18080 

18.59 

II 

0.0087 

173.5 

88.4 

529.1 

58.00 

22290 

22.04 

II 

0.0076 

263.2 

148.1 

527.8 

102.6 

47340 

33.99 

II 

0.0049 

248.3 

138.0 

527.8 

94.67 

42590 

35.9 

II 

0.0052 

233.5 

128.1 

527.8 

89.40 

38110 

33.97 

II 

0.0055 

213.5 

114.8 

527.7 

77.62 

32400 

29.50 

II 

0.0061 

93.4 

101.7 

527.7 

68.13 

27070 

25.89 

II 

0.0068 

277.1 

157.7 

527.7 

109.2 

51930 

41.51 

II 

0.0046 

Ka 

(md . ) 
225.90 
231.95 
222.28 
213.74 

204.51 
176.30 
162.06 

160.51 
142.36 
134.54 
112.62 
102.86 

98.35 

100.82 

106.27 

108.41 

113.80 

95.59 


. 

i.ew 

x.e  s 

Run 

No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


-  82  - 


TABLE  15 
GAS  FLOW  DATA 

CORE  #9  INDIANA  LIMESTONE 


Pressure-Psia 

Pi  P2 

Tflow  Qc 

(average) 

°R  (SCFH) 

(P!2-P22) 

(Psia)2 

G 

(lbs  ft 

2sec) 

1 

Pavg 

(Psi-1) 

Ka 

(md. 

24.45 

13.72 

532 

0.0374 

410 

0.0142  x 

10‘2 

0.0524 

4.01 

17.79 

13.62 

530 

0.0128 

131 

0.0049 

If 

0.0637 

4.25 

24.36 

13.62 

530 

0.0376 

408 

0.0143 

•i 

0.0527 

4.03 

38.82 

13.63 

530 

0.114 

1321 

0.0433 

If 

0.0381 

3.76 

45.6 

13.63 

530.1 

0.151 

1889 

0.0574 

ii 

0.0338 

3.50 

113.85 

13.73 

531.2 

0.891 

12750 

0.339 

ft 

0.0157 

3.08 

156.0 

13.89 

531.2 

1.598 

24140 

0.607 

II 

0.0118 

2.91 

216.0 

14.34 

531.7 

2.87 

46450 

1.091 

II 

0.0087 

2.73 

263.6 

14.57 

532.2 

4.08 

69270 

1.550 

It 

0.0072 

2.62 

310.0 

15.57 

532.2 

5.42 

95860 

2.060 

II 

0.0061 

2.51 

395.5 

16.77 

528.6 

7.02 

128920 

2.668 

II 

0.0053 

2.41 

412.5 

18.47 

529.4 

8.91 

169820 

3.386 

If 

0.0046 

2.32 

459.0 

20.28 

529.4 

10.75 

210270 

4.085 

II 

0.0042 

2.27 

511.7 

22.67 

529.4 

12.91 

261320 

4.906 

II 

0.0037 

2.19 

558.7 

25.15 

529.4 

14.90 

311510 

5.662 

II 

0.0034 

2.12 

609.7 

28.12 

528.7 

17.25 

370940 

6.555 

ft 

0.0031 

2.07 

707.0 

35.56 

529.9 

21.67 

498530 

8.235 

If 

0.0027 

1.95 

670.3 

32.82 

529.6 

20.04 

448160 

7.615 

If 

0.0028 

2.00 

642.9 

31.01 

529.6 

18.76 

412360 

7.129 

II 

0.0030 

2.04 

687.1 

27.44 

529.6 

20.74 

471290 

7.882 

1 1 

0.0028 

1.97 

*0d,0 

Run 

No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 
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TABLE  16 
GAS  FLOW  DATA 
CORE  if/ 11  BASAL  QUARTZ 


Pressure 

pi 

-Psia 

P2 

low 

(average) 

°R 

(SCFH) 

2  2 
<P1  p2  > 

(Psia)2 

G 

(lbs/ft2sec 

20.50 

13.51 

'  536.4 

2.53 

238 

0.924  x  10 

29.97 

13.89 

536.4 

7.20 

705 

2.630 

40.76 

15.12 

535.0 

13.95 

1433 

5.096 

55.0 

17.83 

534.9 

22.98 

2709 

8.395 

74.9 

24.71 

534.9 

40.05 

5002 

14.63 

95.0 

34.03 

534.6 

56.13 

7871 

20.51 

113.0 

43.2 

534.6 

71.87 

10900 

26.26 

133.3 

53.7 

534.6 

89.54 

14880 

32.71 

213.7 

97.7 

534.4 

163.4 

36120 

59.69 

192.5 

86.0 

534.4 

144.50 

29640 

52.79 

173.2 

75.3 

534.4 

125.9 

24330 

45.98 

152.0 

64.3 

534.4 

108.1 

18960 

39.49 

' 

■ 

Run 

No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 
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TABLE  17 
GAS  FLOW  DATA 
CORE  #10  BASAL  QUARTZ 


Pressure 

pi 

-Psia 

P2 

■^f  low 
(average) 
°R 

Qc 

(SCFH) 

(P^-Pa2) 

(Psia)2 

G 

(lbs /ft2 

sec) 

20.28 

13.7 

535.8 

1.347 

225 

0.492  x 

10’ 

30,69 

13.80 

535.8 

4.38 

751 

1.601 

tt 

40,83 

14.2 

534.6 

8.29 

1466 

3,031 

ii 

63,8 

16.4 

534.6 

18.46 

3800 

6.748 

ii 

83,6 

20.3 

534.6 

29.91 

6585 

10,930 

n 

103.9 

26.3 

534.6 

42 ,36 

10110 

15.49 

ii 

125.1 

33.9 

534.6 

55.37 

14510 

20.24 

ii 

143.8 

41.1 

534.2 

68.10 

18990 

24,90 

ii 

166.2 

50.0 

533.9 

83.07 

25130 

30.37 

ti 

183.4 

59.0 

533.8 

94.65 

30150 

34.60 

ii 

203.4 

65.2 

533.1 

110.40 

37130 

40.37 

M 

225.3 

74,4 

532.0 

125,30 

45240 

45.81 

ii 

243.7 

82,4 

531.7 

138.30 

52590 

50.57 

If 

262.6 

90.6 

531.0 

151.8 

60750 

55.44 

ii 

. 


■ 


•.  :  i 

• 

' 

s’  i’ 

-•.e.  ■  ■ 

. 

X  * 

..  : 

s 

.  ' r 

•v 

a.1  e 

■> 

Run 

No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 
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TABLE  18 
GAS  FLOW  DATA 


CORE  #12 

VIKING 

SANDSTONE 

Pressure 

-Psia 

^  f  1  ow 
(average) 

Qc 

2  2 
(pi2*p2  ) 

G 

r\ 

?! 

P2 

°R 

(SCFH) 

(Psia) 

(lbs /ft 

^sec) 

80.2 

13.61 

533.2 

0.04 

6238 

1.462 

x  10" 

148.5 

13.61 

533.6 

0.126 

21860 

4.605 

m 

228.1 

13.61 

534.4 

0.272 

51830 

9.94 

ii 

322.1 

13.60 

534.4 

0.491 

103560 

17.94 

n 

462.9 

13.60 

535.2 

0.917 

214060 

33.51 

ti 

601.8 

13.61 

535.4 

1.424 

361930 

52.04 

ii 

754.5 

13.65 

534.1 

2.051 

569100 

74.95 

II 

897.8 

13.68 

535.8 

2.706 

805900 

98.89 

ii 

1050.7 

13.72 

535.9 

3.460 

1103750 

126.4 

II 

1201.5 

13.77 

535.4 

4.225 

1443400 

154.4 

II 

1298.9 

13.81 

535.4 

4.722 

1686900 

172.6 

ii 

1384.3 

13.85 

534.75 

5.194 

1916100 

189.8 

II 

1481.3 

13.91 

534.75 

5.721 

2194000 

209.1 

ft 

Run 

No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 
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TABLE  19 
GAS  FLOW  DATA 
CORE  #13  VIKING  SANDSTONE 


Pressure-Psia 

P1  P2 

flow 

(average) 

°R 

(SCFH) 

2  2 
(pi  -P2  > 

(Psia)2 

G 

(lbs  'ft2 

67.5 

13.60 

533.7 

0.105 

4364.2 

3.835  x 

116.7 

13.60 

533.7 

0.283 

13430 

10.34 

262.3 

13.59 

532.85 

1.082 

68620 

39.52 

367.7 

13.62 

533.1 

1.853 

135000 

67.69 

508.5 

13.67 

533.0 

3.081 

258400 

112.6 

664.4 

13.77 

533.0 

4.613 

441200 

168.5 

818.2 

13.93 

531.65 

6.254 

669300 

228.4 

879.6 

14.01 

531.3 

6.934 

773400 

253.2 

1009.6 

14.21 

530.9 

8.477 

1019000 

309.7 

1119.8 

14.36 

529.4 

9.86 

1253800 

360.2 

1217.0 

14.59 

527.8 

11.21 

1481000 

409.5 

1308.0 

14.83 

527.5 

12.35 

1711000 

451.1 

1409.2 

15.16 

527.15 

13.71 

1986000 

500.9 

1493.6 

15.46 

527.15 

14.78 

2231000 

539.9 

. 

»»• 

Run 

No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 
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TABLE  20 
GAS  FLOW  DATA 

CORE  #14  BELLY  RIVER  SANDSTONE 


Pressure-Psia 

pl  P2 

Tflow 

(average) 

°R 

(SCFH) 

(,t]2-?22'> 

(Psia)^ 

G  9 

(lbs/f tzsec) 

1 

Pavg 

(Psi"1) 

Ka 

(md . ) 

1428.9 

16.84 

530.4 

19.1 

2041600 

6.980  x 

io"2 

.0014 

0.430 

802.2 

14.14 

530.6 

7.298 

643300 

2.668 

m 

.0025 

0.510 

899.6 

14.33 

530.6 

8.849 

809100 

3.235 

II 

.0022 

0.493 

1006.6 

14.65 

530.6 

10.68 

1013000 

3.9050 

n 

.0020 

0.477 

1098.1 

14.99 

529.6 

12.32 

1206000 

4.505 

II 

.0018 

0.462 

424.5 

13.72 

530.3 

2.371 

180000 

.8668 

It 

.00456 

0.584 

514.8 

13.77 

530.8 

3.349 

265000 

1.224 

ii 

.00378 

0.562 

612.2 

13.86 

530.9 

4.538 

374500 

1.659 

II 

.00319 

0.541 

703.9 

13.97 

530.9 

5.788 

495300 

2.116 

II 

.00279 

0.523 

1192.6 

15.43 

530.3 

13.95 

1422000 

5.099 

ii 

.00163 

0.446 

1323.0 

16.11 

529.4 

16.460 

1750000 

6.017 

•t 

.00149 

0.429 

309.1 

13.71 

530.9 

1.328 

95380 

0.4855 

II 

.00620 

0.615 

119.5 

13.68 

532.7 

0.236 

14090 

0.0863 

II 

.0150 

0.738 

92.4 

13.68 

533.0 

0.145 

8350 

0.05297 

II 

.0189 

0.764 

70.1 

13.68 

533.0 

0.0864 

4730 

0.03157 

If 

.0239 

0.805 

44.1 

13.68 

533.7 

0.0344 

1760 

0.01258 

II 

.0346 

0.862 

31.40 

13.68 

533.7 

0.0173 

799 

0.006324 

II 

.0444 

1.0334 

> 

... 

. 
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TABLE  21 

EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 
OF  CONSOLIDATED  POROUS  SAMPLES 

INDIANA  LIMESTONE  LONG  CORE 


Run 

(pl2-  P22)Mgc 

G 

f 

Re 

No . 

LB2ZRT  jd  G 

ILL 

1 

13 

2-195  x  10 

4-64  x  102 

16.9 

0.135 

2 

2- 114 

3-32 

22.7 

0.0965 

3 

2.050  " 

2.27 

32-2 

0  066 

4 

2-005  " 

1.37 

52-3 

0.0398 

5 

1.977 

0.98 

72.0 

0.0285 

6 

1  944 

0.65 

106.8 

0.0189 

7 

1.869  " 

0.40 

167 

0.0116 

8 

1.805 

0.20  " 

322 

0.00582 

9 

1.319 

0  10 

471 

0.00291 

10 

1.948 

0.43 

162 

0.0125 

11 

2.017 

1.12 

64-3 

0.0326 

12 

1.991 

1  00 

71.1 

0.0291 

13 

2.067 

1  -  75 

42-2 

0.0509 

• 

TABLE  22(A) 


EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 
OF  CONSOLIDATED  POROUS  SAMPLES 

(BEFORE  PRESSURE  TAPS  DRILLED) 

ALUNDUM  CORE  (TA7921) 


(P12-P22)Mgc 

LB2ZRT^OG 

G 

~F 

£ 

Re 

1.435  x  1011 

534.3 

20.9 

0.0985 

1.572 

923.1 

13.2 

0.170 

1.671 

1223.4 

10.6 

0.225 

1.781 

1482.4 

9.33 

0.273 

1.836  « 

1693.5 

8.42 

0.312 
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TABLE  22(B) 

EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 
OF  CONSOLIDATED  POROUS  SAMPLES 


ALUNDUM  CORE  (TA  7921) 

(PRESSURE  TAPS  DRILLED  ALONG  LENGTH  OF  CORE. 
PRESSURE  TAPS  UPSTREAM  AND  DOWNSTREAM  OF  END  PLATES.) 


Run 

(p^-p,2 

)  Mgc 

G 

f 

Re 

No. 

TZ  2zrT 

D 

^G 

1 

2.278  x 

10U 

66.766  x 

102 

4.64 

.756 

2 

2.115 

tt 

56.395 

it 

5.10 

.639 

3 

1.763 

It 

32.069 

ft 

7.48 

.363 

4 

1.947 

ft 

44.804 

it 

5.91 

.508 

5 

1.544 

tt 

15.822 

it 

13,3 

.179 

6 

1.518 

It 

9.805 

tt 

21.1 

.111 

7 

1.549 

It 

5.134 

it 

41.1 

.0582 

8 

1.497 

It 

2.548 

tt 

80.0 

.0289 

9 

1.180 

M 

1.371 

tt 

117 

.0155 

10 

1.185 

it 

1.353 

it 

119 

.0153 

11 

1.187 

tt 

1.241 

tt 

130 

.0141 

12 

1.192 

it 

0.873 

it 

186 

.00989 

13 

1.018 

ft 

0.043 

tt 

322 

.00487 

14 

1.227 

tt 

3.096 

it 

539 

.00351 

15 

1.169 

tt 

1.732 

it 

919 

.00196 

16 

1.202 

ft 

0.462 

it 

354 

.00052 

17 

1.277 

tt 

5.170 

tt 

33.6 

.0586 

18 

1.257 

it 

4.060 

it 

42.1 

.046 

19 

1.257 

it 

1.970 

tt 

86.8 

.0223 

20 

1.555 

it 

1.808 

it 

117.1 

.0205 

21 

1.457 

tt 

1.253 

it 

158 

.0142 

22 

1.410 

tt 

6.994 

it 

27.4 

.0792 

23 

1.079 

it 

1.704 

tt 

86.2 

.0193 

24 

1.423 

tt 

5.138 

it 

37.7 

.0582 

25 

1.478 

tt 

12.654 

tt 

15.9 

.143 

26 

1.689 

tt 

27.128 

It 

8.47 

.307 

27 

1.736 

ft 

30.375 

ft 

7.78 

.344 

28 

2.015 

tt 

48.186 

ft 

5.69 

.546 

29 

1.887 

it 

40.435 

It 

6.78 

.458 

30 

1.710 

it 

27.613 

If 

8.43 

.313 

31 

1.483 

ft 

8.124 

it 

24.8 

.0920 

■■  ■ 


r.e.i 

II 

II 
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TABLE  22 (Cl) 

EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 
OF  CONSOLIDATED  POROUS  SAMPLES 


ALUNDUM  CORE  (TA  7921) 
(PRESSURE  TAPS  IN  END  PLATES) 


Run 

(Pi2-p22 

)  Mgc 

G 

f 

Re 

No 

Lb  2ZRT 

flG 

h 

1 

1.293  x 

1011 

1.376  x 

10J 

564 

.0036 

2 

1.313 

1! 

1.834 

ti 

430 

.0048 

3 

1.286 

II 

2.497 

it 

309 

.00653 

4 

1.272 

II 

1.898 

it 

403 

. 00496 

5 

1.191 

It 

0.110 

it 

6503 

.00029 

6 

1.196 

It 

0  165 

ii 

4353 

.000431 

7 

1.226 

It 

0.208 

ii 

3540 

. 000544 

8 

1.188 

It 

0.068 

it 

10490 

.000178 

9 

1.199 

ft 

0.394 

tt 

1828 

.00103 

10 

1.278 

II 

0.609 

If 

1260 

.00159 

11 

1.286 

It 

0.917 

ii 

842 

.0024 

12 

1.190 

II 

0.0154 

If 

46410 

.0000403 

13 

1.177 

II 

0.0079 

n 

89480 

.0000217 

14 

1.016 

II 

0.0026 

ii 

234700 

.0000068 

15 

1.273 

II 

1.844 

VI 

415 

.00482 

16 

1.282 

II 

2.446 

It 

315 

.00639 

17 

1.308 

If 

2.897 

II 

271 

.00757 

18 

1.288 

it 

4.179 

If 

185 

.0109 

19 

1.279 

it 

5.329 

II 

144 

.0139 

20 

1.284 

It 

6.284 

It 

123 

.0164 

21 

1.288 

it 

7.306 

II 

106 

.0191 

22 

1.321 

ft 

8.287 

II 

95.1 

.0217 

23 

1.293 

it 

9.061 

II 

85.7 

.0237 

24 

1.303 

tt 

9.498 

ft 

82.4 

.0248 

25 

1.287 

ii 

7.900 

II 

97.8 

.0207 

26 

1.313 

tt 

10.903 

II 

72.3 

.0285 

27 

1.271 

IV 

5.365 

If 

127 

.014 

' 


-v* 
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ti 

.  . 

t'Oi.O 

it 

' 

' 

ce.es 

* 

eeaoo . 

VuiC’-. 

Run 

No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 
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TABLE  22 (C2) 

EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 
OF  CONSOLIDATED  POROUS  SAMPLES 

(PRESSURE  TAPS  IN  END  PLATES) 

ALUNDUM  CORE  (TA7921) 


(Pl2-P,2)Mgc 

G 

Lg2ZRT  jA*  G 

P* 

1.191  x  1011 

2.766  x  10 

1.193 

4.561  " 

1.205 

6.566  " 

1.273  " 

14.80  " 

1.264 

27.84  " 

1.290 

35.94  " 

1.285 

46.93 

1.314 

59.16  " 

1.326 

73.24 

1.368 

87.25 

1.364 

79.39 

1.369 

64.83 

Re 


(tSf^AT)  HJ100  MUdPfUJA 


" 

• 

.  . 

Run 

No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 
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TABLE  23 


SUMMARIZED  CALCULATIONS  OF  THE 
EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 
OF  CONSOLIDATED  POROUS  SAMPLES 


CORE  #1 

(Pl2~P?2)  Mgc 

Lb  2ZRT  f^G 

CARDIUM 

G 

SANDSTONE 

f 

Re 

4.053  x  1012 

6.767 

x  102 

66.4 

.0284 

4.192 

9.816 

it 

47.3 

.0412 

4.217 

13.32 

If 

35.1 

.0559 

4.381 

19.65 

it 

24.7 

.0825 

4.536 

26.97 

it 

18.6 

.113 

4.658 

32.36 

it 

15.9 

.136 

4.861 

39.89 

it 

13.5 

.167 

5.007  " 

45.64 

ii 

12.2 

.192 

5.178 

51.69 

n 

11.1 

.217 

5.335 

56.74 

n 

10.4 

.238 

5.514 

62.81 

19 

9.73 

.264 

5.612 

66 . 36 

II 

9.37 

.279 

5.308 

66.50 

it 

8.84 

.279 

5.667 

83.14 

ti 

7.55 

.349 

6.027 

101.01 

it 

6.61 

.424 

6.541 

125.25 

ii 

5.79 

.526 

7.023  " 

151.15 

ii 

5.15 

.635 

7.488 

174.94 

it 

4.74 

.734 

7.883 

205.34 

ii 

4.25 

.862 

10.779 

360.09 

ii 

3.32 

1.515 

10.255 

329.59 

ii 

3.45 

1.384 

9.706 

302.55 

ii 

3.55 

1.270 

9.178  " 

270.48 

ii 

3.76 

1.135 

8.875  " 

237.84 

ti 

4.13 

0.998 

8.021 

212.11 

II 

4.19 

.890 

■  . 


■  1 

> 

. 

• 
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TABLE  24 

EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 

OF  CONSOLIDATED  POROUS  SAMPLES 

CORE  #2  CARDIUM  SANDSTONE 


Run 

No 

(px2  p,2 

Tg  2ZET 

)  Mgc 

jl/G 

G 

f 

Re 

1 

0.695  x 

1016 

0.993 

130.8 

00552 

2 

1.040 

ii 

1.106 

176.6 

.00616 

3 

1.119 

ii 

1.550 

134.9 

.00863 

4 

1.470 

ii 

2.450 

112.1 

.0136 

5 

1.624 

ti 

4.053 

74.9 

.0226 

6 

1.694 

ii 

3.828 

82,7 

.0213 

7 

1.831 

If 

6.093 

56.2 

.0339 

8 

1.939 

ti 

9.021 

40.2 

.0502 

9 

2.017 

ii 

12.119 

31.1 

.0675 

10 

2.088 

n 

16.873 

23.1 

.0939 

11 

2.174 

ii 

23.176 

17.5 

.129 

12 

2.238 

II 

29.244 

14.3 

.163 

13 

2.303 

it 

35.398 

12.2 

.197 

14 

2.399 

it 

44.110 

10.2 

.246 

15 

2.496 

n 

52.452 

8.89 

.292 

16 

2.485 

II 

59.294 

7.83 

.330 

17 

2.637 

II 

66.154 

7.45 

.368 

18 

2.717 

it 

75.709 

6.71 

.421 

:  \{  )  V  ,  I'-' 

•  r  .  '  ; 


n  >  . 

•I 

c  £ 

di 
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TABLE  25 

EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 
OF  CONSOLIDATED  POROUS  SAMPLES 


CORE  #3  CARDIUM  SANDSTONE 


Run 

No 

(Pl2-P,2)Mgc 

'  Lb  zzki^g 

G 

f 

Re 

1 

4.785  x  1012 

0.0032  x 

104 

1181 

.00143 

2 

4.968  " 

0.0107 

11 

366 . 6 

.00478 

3 

5.23 

0.0324 

11 

127.5 

.0145 

4 

5.29 

0.0764 

11 

54.7 

.0342 

5 

6.03 

0.145 

ii 

32.8 

.0648 

6 

6.06 

0.275 

11 

17.4 

.123 

7 

6.36 

0.347 

ii 

14.5 

.155 

8 

6.99 

0.510 

IV 

10.8 

.228 

9 

7.76 

0.659 

11 

9.19 

.295 

10 

8.16 

0.880 

11 

7.32 

.393 

11 

9.13 

1.174 

11 

6.14 

.525 

12 

10.08 

1.486 

11 

5.36 

.  664 

13 

11.18 

1.852 

11 

4.77 

.828 

14 

11.95 

2,138 

ii 

4.41 

.956 

15 

18.50 

4.856 

It 

3.01 

2.171 

16 

17.50  " 

4.292 

ii 

3.22 

1.919 

17 

16.10 

3.727 

11 

3.41 

1 . 666 

18 

14.62 

3.044 

11 

3.79 

1.361 

19 

13.33 

2.607 

11 

4.04 

1.166 

)  VK  (  i 


n 

l  AU  i  «’ 

. 
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TABLE  26(A) 

EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 
OF  CONSOLIDATED  POROUS  SAMPLES 


CORE  #4  CARDIUM  SANDSTONE  (Tf  -  73°F) 


Run 

No 

(P,2-P,2)  Mgc 

Tg  zzRrjyte 

G 

f 

Re 

I 

9.502  x  1012 

4.407  x 

104 

2.32 

2.564 

2 

8.459 

4.013 

ii 

2.27 

2.335 

3 

7.869  " 

3.555 

n 

2.39 

2.068 

4 

7.371 

3.108 

ii 

2.56 

1.808 

5 

6.505 

2.449 

ft 

2.86 

1.425 

6 

5.922 

2.003 

ii 

3.19 

1.165 

7 

5.396 

1.635 

it 

3.56 

0.951 

8 

4.917 

1.298 

ii 

4.08 

0.755 

9 

4.469  " 

0.979 

it 

4.92 

0.570 

10 

4.069  " 

0.664 

ff 

6.61 

0.386 

11 

3.506 

0.393 

it 

9.62 

0.229 

12 

3.203 

0.205 

ii 

16.84 

0.119 

13 

3.072  " 

0.112 

ii 

29.57 

0.0652 

14 

2.246 

0.051 

ft 

47.47 

0.0297 

ntf  ' 

t}l,  £ 

•Vi.  i 

"  eoe.t) 

- 

r, 

.  ■ 

Run 

No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 
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TABLE  26(B) 

EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 
OF  CONSOLIDATED  POROUS  SAMPLES 


CORE  #4  CARDIUM  SANDSTONE 

(Tf 

=  95°F) 

(P,2-P22)Mgc 

G 

f 

Re 

tg  2ZRT^tG 

r- 

3.049  x  1012 

0,141  x 

104 

23.31 

.0820 

2.928 

0.062 

II 

50.91 

.0361 

3.408 

0.331 

II 

11.10 

.193 

4.042 

0.648 

II 

6.72 

.377 

4.716 

1.122 

(I 

4.53 

.653 

3.238 

1.471 

I! 

3.84 

.856 

5.682 

1.756 

II 

3.49 

1.022 

6.018 

2.015 

II 

3.22 

1.172 

6.417 

2.314 

II 

2.99 

1.346 

7.178 

2.877 

II 

2.69 

1.674 

8.835  " 

4.193 

If 

2.27 

2.439 

8.075 

3.611 

II 

2.41 

2.101 

7.587  " 

3.217 

II 

2.54 

1.872 

. 

TABLE  26(C) 


EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 
OF  CONSOLIDATED  POROUS  SAMPLES 


CORE  #4  CARDIUM  SANDSTONE  (Tf  =  115 °F) 


Run 

No 

1 

2 

3 

4 

5 

6 
7 


(P12-P22)  Mgc 

G 

f 

Lg 

P 

4.344  x  1012 

0.830  x  104 

5 . 64 

7.975 

3.539 

2.43 

7.376 

3.025  " 

2.63 

6.783 

2.560 

2.86 

6.169 

2.096  " 

3.17 

5.562 

1.664 

3.60 

4.862 

1.177 

4.45 

4.109 

0.7174  " 

6.26 

Re 

.483 

2.059 

1.760 

1.489 

1.219 

0.956 

0.685 

0.417 


8 


*dd  x 

1 

2 

3 

4 

5 

6 

1 

8 

9 

10 

11 

12 

13 

14 

15 

16 
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TABLE  27 

EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 
OF  CONSOLIDATED  POROUS  SAMPLES 
CORE  #5  CARDIUM  SANDSTONE 


(Pl2-P22)M8c 

G 

f 

Re 

Lfi  ZZRTp-G 

H' 

14.72  x  1012 

4.348  x 

104 

3.13 

2.170 

13-50 

3.805 

Tf 

3.28 

1.899 

15.73 

4.808 

IV 

3.02 

2.399 

15.55 

4.783 

tl 

3.00 

2.387 

12.63 

3.321 

II 

3.51 

1.657 

11.20 

2.735 

II 

3.78 

1.365 

10.45 

2.328 

ti 

4.15 

1.162 

9.47 

1.961 

II 

4 . 46 

.979 

8.61 

1.601 

II 

4.97 

.799 

7.41 

1.128 

II 

6.07 

.563 

6.79 

0.892 

II 

7.03 

.445 

4.75 

0.241 

ti 

18.21 

.120 

5.61 

0.522 

It 

9.93 

.275 

5.02 

0.120 

II 

38.65 

.060 

4.24 

0.0535 

It 

73.21 

.0267 

3.61 

0.0082 

II 

406.7 

.0041 

■  *;'•  i 


• 

.  >.a 
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TABLE  2  8 

EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 
OF  CONSOLIDATED  POROUS  SAMPLES 

CORE  #6  STURGEON  LAKE  SOUTH  LIMESTONE 


Run 

2 

O’!  -*2 

2 

)  Mgc 

G 

f 

Re 

No 

Lg  2ZRT G 

1 

7.270 

KlO13 

0.012  x 

i  ^ 

o 

r—t 

4637 

•00017 

2 

7.838 

f  1 

0.027 

tl 

2222 

.00038 

3 

8.604 

u 

0.056 

II 

1176 

.00079 

4 

9.246 

ti 

0.087 

11 

813 

.00123 

5 

10.152 

n 

0.133 

II 

584 

.00187 

6 

10.448 

it 

0.173 

II 

462 

.00244 

7 

10.805 

Tl 

0.230 

f  1 

360 

.00324 

8 

11.257 

11 

0.318 

II 

271 

.00448 

9 

11.916 

tl 

0.422 

It 

216 

.00595 

10 

11.983 

II 

0.508 

II 

181 

.00716 

11 

16.153 

II 

3.490 

tl 

35.4 

.00492 

12 

17 .025 

II 

5.753 

II 

22.7 

.00811 

13 

18.467 

11 

7.064 

11 

20.0 

.00995 

14 

19.304 

II 

8.720 

11 

16.9 

.0123 

15 

19.840 

II 

10.092 

II 

15  .0 

.0142 

16 

20.240 

II 

10.855 

If 

14.3 

.0153 

17 

20.956 

II 

12.787 

II 

12.5 

.0180 

18 

21.350 

II 

14.154 

11 

11.5 

.0199 

19 

22-200 

11 

16.217 

II 

10.5 

.0228 

20 

29.097 

II 

41.620 

If 

5.35 

.0586 

21 

28.340 

II 

38.025 

it 

5.70 

.0536 

22 

27  .534 

II 

34.583 

11 

6.10 

.0487 

23 

26.913 

II 

31.286 

II 

6.58 

.0441 

24 

25.979 

II 

28.71 

II 

6.93 

.0405 

25 

25.030 

M 

24.12 

II 

7.94 

.0340 

26 

24.064 

II 

20.89 

It 

8.82 

.0294 

27 

23.261 

II 

17.82 

It 

9.99 

.0251 
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TABLE  29 

EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 
OF  CONSOLIDATED  POROUS  SAMPLES 

CORE  #7  STURGEON  LAKE  SOUTH  LIMESTONE 


(Pl2-P22)  Mgc 

LB  2ZRT/*G 

G 

f 

Re 

3.476  x  1013 

5.95 

1075 

.000091 

3.788 

23.64 

295 

.000361 

4.132 

54.7 

139 

.000835 

4.461 

118.6 

69.2 

.00181 

5.680 

468.9 

22.3 

.00716 

6.693 

878.5 

14.0 

.0134 

8.018 

1677.9 

8.79 

.0256 

9.157 

2412.0 

6.99 

.0368 

10.46 

3318.6 

5.80 

.0506 

11.79 

4193.3 

5.17 

.0640 

12.89 

5175-1 

4.58 

.07  90 

14.02 

6064.3 

4.25 

.9025 

15.51 

8135.4 

3.51 

.124 

15.08 

7537.1 

3.68 

.115 

14.63 

6807.0 

3.95 

.104 

' 
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TABLE  30 

EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 
OF  CONSOLIDATED  POROUS  SAMPLES 

CORE  #8  BEREA  SANDSTONE 


2  2 
(p,  -v 

lbTzrT 

)Mgc 

[LG 

G 

f 

Re 

4.134  x 

10U 

0.428  x 

102 

1688 

.00085 

4.026 

fl 

0.205 

19 

3433 

.00041 

4.199 

99 

2.684 

19 

273 

.00534 

4-368 

ft 

10.49 

If 

72-8 

.0209 

4.566 

it 

19.64 

99 

40.6 

.0391 

5.299 

99 

44-05 

ii 

21.0 

.0877 

5.767 

ft 

67.46 

19 

14.9 

.134 

5.823 

tt 

67.005 

91 

15.2 

.133 

6.567 

II 

107.02 

ii 

10.7 

.213 

6.950 

II 

130.89 

19 

9.28 

.261 

7.628 

f  1 

156.54 

99 

8.52 

•  312 

7.928 

ff 

185.34 

19 

cc 

.369 

9.520 

If 

326.82 

99 

5-09 

.651 

9.286 

f  1 

301.76 

99 

5-38 

.601 

8.808 

9  9 

285.28 

19 

5.40 

.570 

8.633 

91 

247.90 

99 

6.09 

.494 

8.223 

9  9 

217  .75 

99 

6.60 

.434 
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TABLE  31 

EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 
OF  CONSOLIDATED  POROUS  SAMPLES 

CORE  #9  INDIANA  LIMESTONE 


(Pl2-P?2)Mgc 

lb  2ZRT  fJLG 

G 

f 

Re 

2.329  x  1013 

0.120  x 

102 

1968 

.000683 

2.191 

0.0412 

M 

5392 

.000235 

2.323  " 

0.1211 

II 

1945 

.00069 

2.480 

0.367 

II 

685 

.00209 

2.670 

0.486 

II 

557 

.00277 

3.044 

2.856 

II 

108 

.0163 

3.196 

5.112 

II 

63-4 

.0291 

3.424 

9.156 

•  I 

37.9 

.0521 

3.581 

12-983 

II 

28.0 

.0739 

3-725 

17.218 

II 

21.9 

.098 

3.890  " 

22.265 

II 

17.7 

.127 

4.034 

28.273 

II 

14.5 

.161 

4.133 

34.044 

II 

12.3 

.194 

4-269 

40.794 

II 

10.6 

•  232 

4.403 

47.003 

II 

9.50 

.268 

4.529 

54-326 

ii 

8.45 

.309 

4.814 

67.931 

II 

7.19 

.387 

4.689 

62.919 

II 

7.56 

.358 

4.605 

58.835 

II 

7.94 

•  335 

4.762 

65.100 

II 

7.42 

.371 
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TABLE  32 

EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 
OF  CONSOLIDATED  POROUS  SAMPLES 

CORE  #10  BASAL  QUARTZ 


(P12-P22)  Mgc 

LB  ZZRT^tG 

G 

f 

Re 

3.600  x  1011 

0.0414  x 

104 

157 

.0119 

3-692 

0.125 

It 

53-4 

.0359 

3.816 

0.254 

II 

27.1 

.073 

4.445 

0.567 

II 

14-2 

.163 

4.752 

0.917 

II 

9.36 

.264 

5.150 

1.298 

II 

7.17 

.373 

5.650 

1.696 

II 

6.02 

.488 

6.016 

2.085 

II 

5.21 

.599 

6.522 

2.540 

II 

4.64 

.730 

6.865 

2.892 

it 

4.29 

.831 

7.253 

3.373 

II 

3.89 

.970 

7.816 

3.833 

it 

3.68 

1.102 

8.218 

4.221 

II 

3.52 

1.214 

8.674 

4.629 

II 

3-39 

1.331 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 
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TABLE  33 

EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 
OF  CONSOLIDATED  POROUS  SAMPLES 

CORE  #11  BASAL  QUARTZ 


(P12-P22)  Mgc 

“lb  TzWJjOcr 

G 

f 

Re 

2.184  x  10U 

0.078  x  104 

98.4 

.0196 

2-275 

0.221 

36.2 

.0556 

2  393 

0.428 

19.7 

.108 

2.744 

0.705 

13.7 

.177 

2.906 

1.228 

8.32 

.309 

3-261 

1.719 

6.67 

.433 

3.524 

2.198 

5  •  64 

.553 

3.858 

2.736 

4.96 

.689 

5-113 

4-969 

3-62 

1.251 

4-754 

4-404 

3.80 

1.108 

4-480 

3-838 

4.10 

.966 

4.069 

3-300 

4.34 

.831 

i 1 

• 
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TABLE  34 

EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 
OF  CONSOLIDATED  POROUS  SAMPLES 

CORE  #12  VIKING  SANDSTONE 


Run 

No 

(pl2-P22  )Mgc 

Lg  2ZRT 

G 

f 

Re 

1 

3.457  x  1014 

12.3 

286 

.00505 

2 

3-834 

38.6 

101 

.0158 

3 

4.191 

83.0 

51.5 

.0341 

4 

4.625 

149.4 

31.5 

.0613 

5 

5.084 

277.4 

18.7 

.114 

6 

5.505 

428.3 

13.1 

.176 

7 

6.007 

614.4 

9.96 

.252 

8 

6.380 

804.6 

8.08 

.330 

9 

6.803 

1024.2 

6.77 

.421 

10 

7.253 

1243.1 

5.95 

.510 

11 

7.563 

1385  -6 

5-56 

.569 

12 

7.799 

1520.9 

5.23 

.624 

13 

8.077 

1668.9 

4.93 

.685 
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TABLE  35 

EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 
OF  CONSOLIDATED  POROUS  SAMPLES 

CORE  #13  VIKING  SANDSTONE 


(P,2-P?2)Mgc 

G 

f 

Re 

8.883  x  1013 

32-3 

161 

.00673 

10.126 

86.8 

68.3 

.0181 

13.508  " 

330.6 

23.9 

.0689 

15.461 

563.9 

16.0 

.118 

17.705 

932.2 

11.1 

.194 

20.107 

1388.4 

b  •  48 

.289 

22.467 

1872.9 

7.02 

.390 

23.409 

2072.0 

6.61 

.432 

25-156 

2526.1 

5-83 

.526 

26.623 

2930.2 

5.32 

•  611 

27.706 

3326.0 

4.88 

.693 

28.986 

3654.2 

4.64 

.762 

30.208 

4043-2 

4.37 

.843 

31-386 

4343-7 

4.23 

.905 
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TABLE  36 

EXPERIMENTALLY  DETERMINED  FLOW  CHARACTERISTICS 
OF  CONSOLIDATED  POROUS  SAMPLES 

CORE  #14  SALT  &  PEPPER  SANDSTONE 


Pl^-P22)Mgc 

Lg  2zrt(ug 

G 

f 

Re 

14 

2.184  x  10 

5.609  x 

103 

6 . 16 

0.452 

1.840 

2.192 

II 

13.28 

0.177 

1.903 

2.649 

II 

11.36 

0.214 

1.965  " 

3-185 

If 

9.76 

0.257 

2.029 

3.669 

II 

8.75 

0.296 

1.603 

0.721 

II 

35.16 

0.0581 

1.664 

1.016 

II 

25.90 

0.0819 

1.731 

1.372 

II 

19.96 

0.111 

1.790 

1.744 

II 

16.24 

0.141 

2.101 

4.133 

II 

8.04 

0.333 

2.188 

4.859 

fl 

7.12 

0.392 

1.521 

0.406 

II 

59.24 

0.0327 

1.266 

0.0725 

II 

276.2 

0.00585 

1.222 

0.0446 

II 

433.6 

0.00359 

1.161 

0.0266 

If 

691-3 

0.00214 

1.083 

0.0106 

II 

1618 

0.000853 

0.9783 

0.00533 

(1 

2906 

0.000429 
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CM 

l 

ft) 

o> 

LL. 

If 

I 

O 

K 

<s> 

cn 


CM 

CL 


P 

cn 

M 

CM 


-px!0'3(F  eer!) 


FIGURE  20 


FLOW  DATA  PLOT  TO  OBTAIN  AND  /S 
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FIGURE  21 


FLOW  DATA  PLOT  TO  OBTAIN  o<  AND  B 
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FLOW  DATA  PLOT  TO  OBTAIN  <*  AND  /S 
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FIGURE  23  FLOW  DATA  PLOT  TO  OBTAIN  «  AND  /8 
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EQUATIONS 


(Pi  -I-  P2) 


Qm  =  flow  rate  at  mean  pressure 
and  flowing  temperature 
Qc  =  flow  rate  at  standard 

conditions,  14.65  psia  and 
60°F. 

P^  =  pressure  base  =  14.65  psia 

Tjt  =  average  flowing  temperature 
in  °R 

Tb  =  temperature  base  =  60°F. 

Pi  =:  upstream  pressure 
?2  -  downstream  pressure 


IC 

a 


^m 

A 


0.25573 


T.c 

j. 


(p2) 


A 


K  =  darcies 
a 

=  cp 

Q  =  cc/sec 

L  =  cm 

P  =  atm 

*  2 
A  =  cm 

T  =  °R 


Ka  = 
Tf  = 

L 

A 

P 

Qc  = 


millidarcies 

°R 

micropoise 

in 

.  2 
m 

psia 

SCFH 


.  M  %SiioiJxbnoo 
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M 


gc 


2  (ZRT) ^  L 


=  52.0956  x  10 


lb 


m 


ZT  c  jA.li  lb_  sec 


Z 

Tf 

h 

L 


=  compressibility 


_  o 


R 

micropoise 

in 


Q 


A 


A(p2) 


Q 


.04 


£ 


A 


lb 


m 


ft  sec 


=  20736 


2  2 

.A(P^)  lb^  sec 

G  lb  ft2 

m 


Q  = 

P  = 

A  = 
L  = 


SCFH 


lb  /ft 


in 
psia 


G 

h 


ft 


-l 


—  =  148.8  x  10 


r 


r 


lb 


m 


Q, 


=  Q  x  Meter  Correction  factor  x 


M 


f t^  sec 


micropoise 

x  35.495 


T 


M 


P 

M 


P 

Atm. 


Q  =  CFII 

PM  =  Meter  inlet  pressure  in 

psia 

Tj^  =5  Meter  inlet  temperature 

p  -  Partial  Pressure  of  water 

meter  inlet  (psig) 


■ 
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Sample  Calculation 

Measured  Flow  Data 

Core  #13  -  Viking  Sandstone 


Area  =  0.8101  in^  Length  =  1.4669  in  Porosity  =  10.28% 

Volume  Rate  =  2.113  CFH 

Temperature  (meter)  =  533.1  °R  Vapor  pressure  =  0.8203  in  Hg 

Pressure  (atmospheric)  =  73.3  cm  Hg  =  13.567  psi 

Flowing  fluid  -  nitrogen 

Temperature  (average  flowing)  =  533.1  °R 

/*tavg.  T.  and  P.)  =  178’6  “icropoise 
P^  =  367.7  psia 

P2  =  13.62  psia 

Pbase  =  14-65  psia 

T,  =  520  °R 

base 


Calculations 


1#  Q(SCFH)  Q (meter) 


P  T 

meter  base 

x  Meter  Correction  Factor  x  ™  x 


T 


meter  ^base 


2. 


3. 


Q(SCFH)  =  X-853  SCFH 


G  =  0.04  0  ( SCFK)/7  =  1.694  X  10 

A  / 


-2 


lb 


ra 


A(p2) 


A 


=  20736 


ft  sec 


A(P2) 


4.136  x  10' 


11  lb  "  sec 


lb  ft* 
iUm  ^ L 


4. 


Mg 


2(ZRT)yU  L 

5.  Mg  c  A(P2) 
2  (ZRT)  LG 


=  52.0956  x  10  x 


13 


15.461  x  10  feet 


Z  ULTf 
-2 


=  3.738  x  10' 


lb 


m 


ft  sec 
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6.  G  5  G  -1 

77,  ~  148.8  x  10  77  =  564  feet 


7. 


N 


Re 


=  0.118 


8. 


f  = 


A(P  )  Mgc 

* 

^(ZRT)  LG 


=  16.0 


/'■  . 
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APPENDIX  III 


1-02 
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FIGURE  26  COMPRESSIBILITY  FACTOR  OF  NITROGEN 

(from  Sage  and  Lacey  26) 
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figure  27  viscosity  of  nitrogen 

(from  Michels  and  Gibson  ^2  ) 


FIGURE  28  VAPOUR  PRESSURE  OF  WATER 

(Chem-  Eng  Handbook ) 
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FIGURE  29  RELATIVE  DENSITY  OF  MERCURY 

l from  Handbook  of  Chemistry  and  Physics 
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FIGURE  30  RELATIVE  DENSITY  OF  WATER 


